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ABSTRACT 

One hundred fifty-four discrete non-nuclear Ultra-Luminous X-ray (ULX) sources, with 
spectroscopically-determined intrinsic X-ray luminosities > 10 39 ergs s _1 , are identified in 82 galaxies 
observed with Chandra's Advanced CCD Imaging Spectrometer. Source positions, X-ray luminosities, 
and spectral and timing characteristics are tabulated. Statistical comparisons between these X-ray 
properties and those of the weaker discrete sources in the same fields (mainly neutron star and stellar- 
mass black hole binaries) are made. Sources above ^10 38 ergs s _1 display similar spatial, spectral, 
color, and variability distributions. In particular, there is no compelling evidence in the sample for a 
new and distinct class of X-ray object such as the intermediate-mass black holes. 83% of ULX candi- 
dates have spectra that can be described as absorbed power laws with index (r) = 1.74 and column 
density (Nh) — 2.24 x 10 21 cm~ 2 , or ^5 times the average Galactic column. About 20% of the ULXs 
have much steeper indices indicative of a soft, and likely thermal, spectrum. The locations of ULXs in 
their host galaxies are strongly peaked towards their galaxy centers. The deprojected radial distribu- 
tion of the ULX candidates is somewhat steeper than an exponential disk, indistinguishable from that 
of the weaker sources. About 5-15% of ULX candidates are variable during the Chandra observations 
(which average 39.5 ks). Comparison of the cumulative X-ray luminosity functions of the ULXs to 
Chandra Deep Field results suggests ~25% of the sources may be background objects including 14% 
of the ULX candidates in the sample of spiral galaxies and 44% of those in elliptical galaxies implying 
the elliptical galaxy ULX population is severely compromised by background active galactic nuclei. 
Correlations with host galaxy properties confirm the number and total X-ray luminosity of the ULXs 
are associated with recent star formation and with galaxy merging and interactions. The preponder- 
ance of ULXs in star-forming galaxies as well as their similarities to less-luminous sources suggest they 
originate in a young but short-lived population such as the high-mass X-ray binaries with a smaller 
contribution (based on spectral slope) from recent supernovae. The number of ULXs in elliptical 
galaxies scales with host galaxy mass and can be explained most simply as the high-luminosity end 
of the low-mass X-ray binary population. 

Subject headings: galaxies: general — surveys — X-rays: binaries — X-rays: galaxies — X-rays: 
general 



1. INTRODUCTION 

Among the most intriguing objects in the X-ray 
sky are the discrete non-nuclear Ultra-Luminous X-ray 
sources (ULXs) in nearby galaxies. This name describes 
sources considerably more luminous than expected for a 
spherically- accreting object of typical neutron star mass. 
Here, ULXs are defined to be those with apparent (i.e., 
assumed isotropically emitting) intrinsic luminosities in 
excess of 10 39 ergs s _1 in the 0.5-8.0 keV bandpass. 

Einstein X-ray images of galaxies revealed 16 ULX can- 
didates (Long & Van Speybroeck 1983; Helfand 1984; 
Fabbiano 1989). ROSAThas extended this list to nearly 
100 (Roberts & Warwick 2000; Colbert & Ptak 2002). 
Yet these objects are rare, reported at a rate of only 1-2 
per galaxy from pointed X-ray observations (Colbert et 
al. 2004; Kilgard et al. 2002; Colbert & Ptak 2002; Fos- 
chini et al. 2002; Humphrey et al. 2003; Irwin, Athcy, 
& Bregman 2003) and perhaps occurring much less fre- 
quently in the nearby Universe as a whole (Ptak & Col- 
bert 2004). ULXs may represent the high-luminosity end 
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of a continuous distribution of typical X-ray sources such 
as supernovae and X-ray binaries (e.g., Grimm, Glifanov, 
& Sunyaev 2003) or they may include new classes of ob- 
jects including intermediate mass black holes (Colbert & 
Mushotzky 1999; Makishima et al. 2000; Colbert & Ptak 
2002; van der Marel 2004), beamed sources (King et al. 
2001; Georganopoulos et al. 2002; Kording, Falcke, & 
Markoff 2002), and hypernovae (Wang 1999). 

Part of the reason ULXs are, as a class, poorly under- 
stood is that past X-ray observatories lack the combina- 
tion of high angular resolution and moderate spectral res- 
olution needed to adequately characterize discrete X-ray 
sources in nearby galaxies. The highest spatial resolution 
instruments aboard both Einstein and ROSAT provide 
little or no spectral information and their spectroscopic 
instruments have poor angular resolution. The ASCA 
satellite provides moderate spectral resolution but very 
poor angular resolution. As a result, ULXs are difficult to 
resolve from other nearby X-ray sources and from diffuse 
emission in the field; their positions and morphologies 
are poorly known thus impeding follow-up investigations 
in other wavebands; and ULX X-ray spectra and timing 
properties are often sparsely sampled. 

We have undertaken an extensive X-ray spectrophoto- 
metry survey of the discrete source populations in nearby 
galaxy fields anchored upon archival data obtained with 
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the Chandra X-ray Observatory Advanced CCD Imag- 
ing Spectrometer (ACIS). The main goal of the survey 
is to systematically search for and evaluate the physical 
properties of ULXs. The survey encompasses a sample 
of ULXs large enough for meaningful statistical analy- 
ses; the survey is based upon high spatial resolution im- 
ages that provide accurate celestial positions as well as 
medium-resolution broad-band spectral signatures and 
time sampling. The galaxies included in the survey span 
the range of Hubble morphological types and include 
galaxies of various mass, gas content, dynamical state, 
and evolutionary history; thus allowing correlations be- 
tween ULXs and their local environments to be studied. 

The galaxy selection process and the methods of X-ray 
data reduction and analysis are described in The re- 
sulting sample of galaxies, their properties, and the prop- 
erties of the ULX candidates are given in The impli- 
cations of the current survey for contemporary physical 
models for the ULX phenomena is discussed in |2] 

2. METHODS 
2.1. Galaxy Selection 

Chandra/ ACIS, operating in imaging mode, provides 
moderate spectral resolution and throughput at the high- 
est spatial resolution available in an X-ray telescope fa- 
cility. As our primary goals are to obtain accurate ce- 
lestial positions and to perform spectroscopic analysis 
of ULX candidates, ACIS is the ideal instrument for the 
present study. To achieve these goals requires a minimum 
of ~50 X-ray counts per source distributed over several 
spectral energy bins. Thus, our sample includes nearly 
all available non-grating, timed exposure mode, ACIS 
galaxy observations (those listing a galaxy or a super- 
nova as a target in the Chandra approved target lists 4 ) 
publicly-available c. May 2003 with >50 counts expected 
from sources with intrinsic 0.5-8.0 keV luminosities of 
L x =10 39 ergs s _1 . NGC 221, NGC 224, and the Mag- 
ellanic Clouds are the only galaxies purposely excluded 
from the sample. These galaxies are known to contain 
no ULXs and require multiple pointings for full coverage 
(except for the dwarf elliptical NGC 221). There were 
no non-targeted galaxies in the Chandra fields that met 
the selection criterion. 

For each candidate galaxy in the archive, the number 
of counts expected from a ULX was estimated using the 
Portable Interactive Multi-Mission Simulator (PIMMS, 
Mukai 1993); the distance, D, to the galaxy; the du- 
ration, t, of the observation; and assuming an absorbed 
power law spectral shape with spectral index r — 1.6 and 
with a hydrogen column density, Nh , equal to the Galac- 
tic column density along the line of sight. The ACIS 
count rate for this spectral shape is ^10~ 4 c s -1 for a 
source flux of 10~ 15 ergs cm -2 s -1 (for a source imaged 
on the back-side illuminated CCDs in the 0.5-8.0 keV 
band and for a typical Nh/10 20 ~ 1-10 cm -2 . The count 
rate is about 30% lower if the source is imaged on the 
front-illuminated CCDs.) Thus, the requirement of 50 c 
from a source with luminosity Lx = 10 39 ergs s -1 can 
be expressed as i/Z? 2 >60 where D is measured in Mpc. 
Typical ACIS observations are ~50 ks duration; limiting 
the distances to host galaxies to ^30 Mpc (^13 Mpc for 

4 http:/ /asc. harvard.edu/targct_lists 



a 10 ks observation). 

Distances to the host galaxies were obtained from the 
literature. In order of preference, distances are based on 
the Cepheid period-luminosity relation (~5% distance 
uncertainty; see Ferrarese et al. 2000), /-band surface 
brightness fluctuations (SBF, 10% uncertainty), on the 
tip of the red giant branch (10%), or on planetary nebula 
luminosity functions (10%). Distances based on globular 
cluster luminosity functions and K' and i^-band SBF 
methods are somewhat more uncertain (Ferrarese et al. 
2000). In a few cases, distances were based on tertiary 
distance indicators such as the brightest blue stars (see 
de Vaucouleurs 1978; Karachentsev & Tikhonov 1994 for 
definitions) and only recent, well-calibrated, distances 
were deemed reliable. Finally, distances tabulated in the 
catalogue of Tully (1988) were used when no other esti- 
mate was available. This catalogue uses the Virgo infall 
model of Tully & Shaya (1984), a variant of the Tully- 
Fisher relation, that assumes an infall velocity of 300 
km s -1 for the Local Group, H Q = 75 km s -1 Mpc -1 , 
and a Virgo distance of 16.8 Mpc. Distances in the Tully 
catalogue are most reliable for galaxies beyond Virgo. 

The duration of the observations, t, were taken from 
the good-time intervals provided as part of the X-ray 
event lists. Galactic hydrogen column densities along 
the line-of-sight to the host galaxies were taken from the 
HI map of Dickey & Lockman (1990) using the FTOOL 
utility nh available from HEASARC. 

2.2. Basic Galaxy Properties 

Several global properties of the host galaxies were ob- 
tained to aid in data reduction and for subsequent anal- 
ysis. Galaxy morphological type, major isophotal angu- 
lar diameter (=_D25, measured at surface brightness level 
25.0 mag sec -2 in blue light), major-to-minor isophotal 
diameter ratio, position angle of the major axis, galac- 
tic latitude, total (asymptotic) B-band magnitudes cor- 
rected for extinction and redshift, and (B — V) colors are 
all taken from the Third Reference Catalogue of Bright 
Galaxies (de Vaucouleurs et al. 1991, hereafter RC3). 
The absolute B magnitudes, Mb, were derived from the 
apparent magnitudes and the adopted distances and con- 
verted to luminosities, Lb, using the standard photomet- 
ric quantities cited in Zombeck (1990). Galaxy center co- 
ordinates were taken from the NASA/IPAC Extragalac- 
tic Database (NED) and are often based on accurate ra- 
dio measurements. Typically, these values differ by no 
more than a few seconds of arc from those reported in 
RC3. 

2.2.1. Far-Infrared Luminosities 

Far-infrared measurements of the host galaxies made 
by the Infrared Astronomical Satellite (IRAS) were used 
to derive the far-infrared luminosities of galaxies in the 
sample. The total flux between 42.4 and 122.5 /im is ap- 
proximated by 1.26 x 10 -11 (2.58S' 6 o+-S'ioo) ergs cm -2 s -1 
where Seo and Sioo are the total flux densities at the 60 
and 100 /Ltm bands, respectively (Rice et al. 1988). Flux 
densities were taken, in order of preference, from the tab- 
ulations in Ho, Filippenko, & Sargent (1997), Rice et al. 
(1988), Knapp et al. (1989), or from NED. 

2.2.2. Nearest Neighbor Distances 
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The separation between each host galaxy and its near- 
est neighbor was taken from the work of Ho, Filip- 
penko, & Sargent (1997) or directly from NED. We de- 
fine nearest neighbor to be the closest galaxy on the 
sky with a small relative velocity between the pair 
(Av<500 km s _1 ). The projected angular nearest- 
neighbor separation, 9p, quoted in this work is given in 
units of the D25 isophotal angular diameter. 

2.3. X-ray Data Reduction 

For each X-ray observation, events within the D25 
isophote were extracted from Level 2 event lists for anal- 
ysis. For events on CCD S4, the destreak algorithm 5 
was applied to remove charge randomly deposited along 
pixel rows during readout and data from all CCDs were 
cleaned of bad pixels and columns. X-ray sources were 
then located and source spectra and time series for events 
within an ~90% encircled energy radius were extracted 
along with a background spectrum in a surrounding an- 
nulus for all sources with estimated luminosities exceed- 
ing ~10 38 ergs s" 1 . (These estimates were based on the 
number of detected counts and assume the same spectral 
shapes as was used for galaxy selection, H2.ll) . Some de- 
tails of the source finding algorithm are given in Swartz 
et al. (2003). The source and background region im- 
age, spectrum, and light curve were then visually in- 
spected for anomolies such as background flares or other 
sources overlapping the extraction regions and appropri- 
ate adjustments made. Sources located within 5" of the 
host galaxy nucleus are not further considered in this 
work. X-ray source positions were mapped to second- 
generation Digitized Sky Survey blue (roughly 400-550 
nm) images of the galaxies in order to identify and reject 
bright (above ~20 mag) foreground stars based on spa- 
tial coincidence. Further analysis of non- X-ray data is in 
progress and will be reported elsewhere. In the few cases 
where more than one compatible observation of a galaxy 
was available, datasets were merged using FTOOL utili- 
ties to increase signal-to-noise. 

2.3.1. Spectral Analysis 

The XSPEC (v. 11. 2) spectral-fitting package was used 
for analysis of events in PI channels corresponding to 0.5 
to 8.0 keV. Spectra were binned as needed to obtain at 
least 20 counts per fitting bin (before background sub- 
traction) to ensure applicability of the \ 2 statistic. For 
sources resulting in fewer than 5 degrees of freedom, the 
unbinned spectra (without background subtracted) were 
also fit using the C-statistic to check for fit parameter 
consistency. Redistribution matrices and ancillary re- 
sponse files were generated using Chandra X-ray Center 
tools and calibration data (CIAO version 2.3) using the 
most recent gain maps and observation-specific bad pixel 
lists and aspect histories. 

For sources imaged on front-side illuminated ACIS 
CCDs, the algorithm developed by Townsley et al. (2000) 
was applied to partially correct for the effects of charge 
loss and smearing caused by charge-transfer inefficiency. 
The algorithm was applied to the Level 1 event list, a 
new Level 2 file generated (following the CIAO threads), 
and the spectra re-extracted. 

5 http: / /asc. harvard.edu/ciao2. 1 / downloads / scripts / destreak.ps 



The acisabs 6 correction was applied (as a multiplica- 
tive model in XSPEC) in all spectral fits to account for 
the temporal decrease in the low-energy sensitivity of the 
ACIS detectors. Sources with detected count rates in ex- 
cess of 0.1 c frame -1 were treated for effects of pileup 
following the procedures described in Davis (2002) and 
using the event pileup model of Davis (2001). No CTI- 
correction was applied to piled-up spectra. 

For each model, the fit parameters, 0.5-8.0 keV model 
flux, x 2 statistic and number of degrees of freedom, 
and corresponding errors were recorded. All errors are 
extremes on the single interesting parameter 90% con- 
fidence intervals except for derived fluxes (and corre- 
sponding luminosities) where the la errors are quoted. 
XSPEC-determined fluxes (see § 12.3. 2[) . averaged over 
the duration of the observation, were then scaled by the 
inverse of the fraction of the telescope point spread func- 
tion (PSF) within the source extraction region and lu- 
minosities computed using the adopted distances to the 
host galaxies. 

X-ray sources in each galaxy field were sorted by num- 
ber of detected counts. Spectral fits were then made 
beginning with the highest-count sources and continuing 
with lower-count sources until all potential ULX candi- 
dates were examined. This resulted in spectral fits for 
many sources with luminosities <10 39 ergs s . 

In addition to fitting models to the 0.5-8.0 keV spec- 
trum of the highest-count sources, the background- 
subtracted X-ray counts for all extracted sources were 
binned into three broad bands; defined as S (0.5-1.0 
keV), M (1.0-2.0 keV), and H (2.0-8.0 keV); and the 
X-ray colors (M - S)/T and {H - M)/T, where T = 
S + M + H, were constructed following Prestwich et al. 
(2003). 

2.3.2. Determining ULX Status Spectroscopically 

By our definition, a ULX is a discrete source whose 
intrinsic luminosity exceeds 10 ergs s -1 in the 0.5- 
8.0 keV energy band. An absorbed (multiplicative 
acisabs and phabs components) XSPEC-specific power 
law model called pegpwrlw was applied to the X-ray spec- 
tra to obtain the intrinsic luminosity. The normalization 
parameter for this model is the intrinsic flux in a pre- 
scribed bandpass. The uncertainty estimate for this pa- 
rameter then provides the uncertainty in the intrinsic 
flux. The observed flux and corresponding uncertainty 
were obtained using XSPEC's flux command. 

The power law model provides a statistically- 
acceptable fit to most of the spectra examined f ^3.2.3f) . 
For these spectra, the average value of the quantity 
(AF/F)^/N ~ 1.7 ± 0.3; where F is the intrinsic flux, 
AF its uncertainty, and N the number of source counts. 
For spectra poorly fit using an absorbed power law, the 
intrinsic flux was estimated using the flux command ap- 
plied to the best-fitting model with the model absorption 
component temporarily set to zero and the uncertainty 
was estimated as AF = 1.7(F/VN). 

2.3.3. Temporal Analysis 

From the time series of events within the source ex- 
traction regions, the X-ray light curves were grouped into 

6 http:/ /www. astro. psu.edu/users/chartas/xcontdir/xcont. html 
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1000-s bins and % 2 tests performed against the constant 
count rate hypothesis to test for source variability. The 
Kolmogorov-Smirnov (K-S) statistic was also computed 
to test the sources for time variabiliy by comparing the 
cumulative event arrival times, binned at the nominal 
frame time (3.24 s in full- frame mode), to that expected 
for a steady source. 

For the subset of variable sources, as deduced from ei- 
ther of these two tests, power density spectra (PDS) were 
also examined. Using the Leahy normalization (Leahy et 
al. 1983), the average power is 2.0 (as expected from Pois- 
son noise) and fluctuations up to 10-15 are commonly 
seen with a maximum power in one or two frequency 
bins typically between 15 and 20 in the normalized power 
spectra. Fluctuations of this order are due to noise. No 
values above these levels were observed. Therefore, no 
statistically-significant periodicities were obtained from 
analysis of the power spectra. 

2.3.4. Spatial Analysis 

The data were checked to ensure the nominal pointing 
accuracy 7 using the latest alignment files (c. 2002 May 
02). The accuracy of absolute positions in the Chandra 
data have a typical rms radius of 0."6 (Aldcroft et al. 
2000). Refinements to absolute positions conceivably can 
be made by cross-correlating X-ray positions with avail- 
able astrometric catalogues. However, this has proved 
impractical for the following reasons: (1) By design, anal- 
ysis is restricted to the small, source-crowded, areas of 
the sky within the host galaxy's L>25- (2) The sample of 
galaxies are typically far from the Galactic plane where 
field stars are relatively rare. 

Refined estimates of the centroids of the sources were 
made by fitting an elliptical Gaussian to the spatial dis- 
tribution of X-ray events. This method works well for lo- 
cating the centroids of the model Chandra point spread 
functions, at various off-axis angles and (monochromatic) 
energies, available from the Chandra X-ray Center PSF 
library 8 . For ULX candidates (those with of order 100 c 
or higher), we estimate the statistical uncertainty in the 
source positions due to centroiding errors is less than 
0."1. Source celestial positions quoted in this work are 
these centroid-refined positions. 

Radial profiles were extracted from a subset of sources 
and compared to model PSFs available from the Chan- 
dra X-ray Center. The subset of sources was chosen as 
those with >10% of X-ray events within the source ex- 
traction region attributable to background. As the back- 
ground region encircles and is contiguous with the source 
region, sources with high background are either embed- 
ded in diffuse emission (a true background) or are ex- 
tended. Where other sources were nearby, these were 
deleted from the data prior to extracting the radial pro- 
file. No ULX candidate appeared significantly extended 
according to this comparison. 

3. RESULTS 

The resulting sample of galaxies and some of their 
properties (distance, Lb, Lfir, and nearest-neighbor dis- 
tance 9p) are listed in Table 1. Also included is a log of 
the corresponding Chandra observations. In a few cases, 

7 http://asc.harvard.edu/cal/ASPECT/fix_ofTset/fix_ofTsct.cgi 

8 Ttp:/ /cda. harvard.edu/pub/arcftp/caldb/acis_psflib_2. 9. tar 



two Chandra observations of a single galaxy were com- 
bined to increase the signal-to-noise but each observa- 
tion is listed separately in Table 1. The observation log 
includes the fraction, /fov, of the angular size defined 
by the D25 ellipse falling within the instrument field of 
view; an estimate, L m i n , of the luminosity of the weak- 
est source for which a statistically-meaningful spectral fit 
could be made; and citations to previous analysis of the 
X-ray data from the literature. Observations performed 
in sub-array mode are noted in the /fov column of Ta- 
ble 1. The tabulated luminosity, L m i n , corresponds to 
roughly 50 source counts. The source detection limit is 
about 10 c. 

3.1. The Sample of Galaxies 

The distribution of galaxies in the Chandra sample 
over Hubble morphological type is shown in Figure ^ 
There are 18 elliptical galaxies, 9 lenticulars, 46 spirals, 
and 9 irregulars and peculiars for a total of 82 galaxies. 
Typically, in this work, galaxies earlier than SO/a are 
collectively referred to as ellipticals and the remainder 
as spirals. While the full range of morphological types 
is represented in the Chandra sample, the 1:2 ratio of 
cllipticals-to-spirals is higher than the ratio 1:4 in the 
Tully (1988) Nearby Galaxy Catalogue and the ~1:3 ra- 
tio in the RC3. 
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Fig. 1. — Distribution of Hubble morphological types among the 
Chandra sample of galaxies. 



The Chandra sample of galaxies is not expected to be 
complete to any specified limiting distance, magnitude, 
or other property. However, the Chandra sample can 
be compared to, for example, the Tully (1988) Nearby 
Galaxy Catalogue of 2367 galaxies. The distribution 
of absolute blue magnitudes of the Chandra sample of 
galaxies is compared in Figure |2 to those of the 1768 
galaxies in the Tully catalogue with tabulated absolute 
blue magnitudes. The average of the Chandra sample is 
0.8 mag brighter in blue light than is the Tully sample. 
The distributions of galaxy distances for the Chandra 
sample and for the Tully catalogue are shown in Fig- 
ure |3 The average distance of galaxies in the Chan- 
dra sample is 11.5±1.2 Mpc compared to 22.5±0.4 Mpc 
for the Tully galaxies. Similarly, the D25 angular sizes 
are larger for the Chandra sample (average 8. '2 ± 0/8) 
compared to the Tully galaxies (3. '8 ± 0.'5). Thus, the 
Chandra sample is comprised of nearby galaxies of large 
angular size and high blue luminosity compared to the 
Tully catalogue of galaxies. The sample is also expected 
be more luminous in X-ray light as X-ray luminosity cor- 
relates with blue luminosity (e.g., Fabbiano & Trinchieri 
1985; Fabbiano 1989) though the correlation varies with 
morphological type and is nonlinear (Fabbiano & Shap- 
ley 2002). The departure of the Chandra sample from 
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the Tully catalogue is not unexpected as the Chandra 
sample represents an X-ray-selected sample. 




-16 -18 -20 

Absolute Blue Magnitude of Galaxy (mag) 

Fig. 2. — Cumulative distribution of the Chandra sample of 
galaxies (solid curve) compared to galaxies from the Tully (1988) 
catalogue (dotted) over absolute blue luminosity. Also shown 
(dashed) is the cumulative distribution of a subsample of the Chan- 
dra galaxies having the same distribution of absolute blue luminos- 
ity as the Tully catalogue. 
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Distance to Galaxy (Mpc) 

Fig. 3. — Cumulative distribution of the Chandra sample (solid 
curve) and subsample (dashed, see text) of galaxies compared to 
galaxies from the Tully (1988) catalogue (dotted) over distance. 

In order to further assess the extent of any selection 
bias in the Chandra sample, a subsample of galaxies 
having the same distribution of absolute B magnitudes 
as Tully-catalogued galaxies within D=29 Mpc was se- 
lected. In § 13.31 this subset is used along with the full 
Chandra sample of galaxies to study correlations between 
the ULX candidates and properties of their host galax- 
ies. The subsample was chosen by first grouping the 
full sample and the Tully catalogue into 4 bins span- 
ning -22.5<Mb<-14.5 with a 2.0 mag width for each bin. 
(There were too few galaxies in the full Chandra sam- 
ple to afford a smaller binsizc.) The subsample was then 
constructed by randomly selecting galaxies from each bin 
such that the fraction of subsample galaxies in each bin 
equals the fraction of Tully galaxies in that bin. The dis- 
tribution of absolute blue magnitudes and distances of 



the subsample are shown in Figures|2]and|31 respectively. 
The average subsample galaxy distance is 10.2 Mpc and 
the average D25 angular size is 9. '3. The Tully catalogue 
of galaxies within 29 Mpc and with tabulated Mb aver- 
age 16.4 Mpc and 4/2. 

The (B—V) colors of the Chandra sample all fall within 
the range of 0.5 to 1.0 mag. Therefore, although blue 
light is more sensitive to the properties of the stellar 
population of the host galaxy than is visible light, the 
blue luminosity is used here as a proxy for galaxy mass 
with an estimated factor of two uncertainty assuming 
My oc mass. On the other hand, a galaxy's far-infrared 
luminosity is proportional to its recent star formation 
rate (Kennicutt 1998). The distribution of the sample 
of galaxies in blue and far-infrared luminosity is shown 
in Figure 01 Elliptical galaxies tend to cluster near the 
upper left in this figure indicative of a relatively high 
mass but little or no recent star formation. In contrast, 
spiral galaxies span from lower-left to upper-right show- 
ing a general trend of star formation rate proportional to 
galaxy mass. It is also evident that spirals have a higher 
star formation rate per unit mass (blue light) compared 
to the ellipticals. 
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Fig. 4. — Distribution of the Chandra sample of galaxies in blue 
and far-infrared luminosity. 

3.2. The ULX Candidates 

Approximately 3400 X-ray sources were detected above 
a signal-to-noise ratio of 2.8 (resulting in >10 c per source 
and much less than 1 false detection per field) within the 
D25 areas of the 82 galaxies analyzed. Among these are 
^1900 sources with intrinsic luminosities estimated from 
detected counts exceeding ~10 38 ergs s _1 . The source- 
region image, spectrum, and light curve were examined 
for all these sources though only 837 had sufficient counts 
(>50) for further analysis. Detailed spectral fits were 
made to 357 of these sources and 154 were determined 
to have luminosities in excess of 10 39 ergs s _1 in the 0.5- 
8.0 keV energy range. 

We may therefore define three source populations: (1) 
The 837-member group above an estimated 10 38 ergs s 
and with >50 c detected for which X-ray colors have been 
determined and the K-S test applied to the X-ray light 
curves, (2) The 357-member subset for which spectral 
parameters and hence luminosities have been established, 



() 



and (3) the subset of 154 ULX candidates. Table 2 lists 
the celestial positions and X-ray properties of the ULX 
candidates. 

The subset of ULX candidates is complete to the ex- 
tent that the derived luminosities (and hence adopted 
distances) are correct. The completeness level of the 357 
sources with spectral fits made can be estimated as the 
luminosity where their cumulative luminosity function 
begins to flatten. This occurs at Ly^^A x 10 38 ergs s _1 
with about 20% lying below this completeness level. 
These lower-luminosity sources tend to be those from 
nearby galaxies with deep exposures but are otherwise 
not expected to differ from sources above the com- 
pleteness limit. The lowest luminosity of a source 
for which spectral parameters have been established is 
8.9 x 10 37 ergs s _1 . The luminosities for the first group 
were estimated from the number of detected counts in the 
manner described in § 12.11 and are therefore much more 
uncertain than for the other two source populations. The 
completeness level for this population, based on its cu- 
mulative luminosity function is also ~4 x 10 38 ergs s^ 1 . 
More importantly for this population, however, is that 
there are sufficient counts detected for analyzing their X- 
ray colors and timing properties without overwhelming 
statistical uncertainties. Thus, this population is limited 
to sources with >50 c detected. 

3.2.1. X-ray Spectra 

The 3-parameter absorbed power law model provides 
one of the simplest phenomonological descriptions of X- 
ray spectra. This model was applied to the observed 
spectra to constrain the observed and intrinsic luminosi- 
ties of the sources (see also § 12.3.21 above) and to estab- 
lish their basic spectral shapes. Statistically-acceptable 
(>90% confidence) fits were obtained using this model 
for 83% (298 of 357) of sources for which spectral pa- 
rameters have been determined and for 130 (84%) of the 
ULX candidates. This does not imply that a power law 
is the only acceptable model for these spectra nor that 
a power law fully describes the physical mechanism(s) 
responsible for the observed X-ray emission. Indeed, an 
absorbed thermal emission line (mekal) model, for exam- 
ple, results in an acceptable fit to ^78% of the sources 
analyzed. 

The fitted power law indices of the ULX candidates 
with statistically-acceptable power-law fits are compared 
in Figure [5] to the weaker sources. Overall, the distribu- 
tions are very similar and can be described by the same 
Gaussian shape. However, there is a significant popula- 
tion of sources (20 of 130) with r>3 in the ULX popu- 
lation that does not have a counterpart in the sample of 
weaker sources (7 of 168). 

Closer inspection showed that sources with a steep 
power law index often also required a high absorbing col- 
umn density that compensates for an observed spectral 
rollover toward low energies. The average column den- 
sity is 7V///10 2O =38.7 ± 3.0 cm" 2 for ULX candidates 
with F>3 compared to 19.6 ± 1.2 for the T<3 sources. 
The ultimate effect is to promote steep power law sources 
with low observed luminosities into the high intrinsic lu- 
minosity ULX candidate category. In contrast, thermal 
models such as the mekal model have a natural rollover 
at low energies, do not require as large an absorbing col- 
umn density, and result in a lower intrinsic luminosity. 
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Power Law Index 

Fig. 5. — Number of ULX candidates (solid histogram) and non- 
ULX sources for which spectral fits were made (dotted) against 
power law index. The two distributions include only sources for 
which an absorbed power law model provides an acceptable fit. 
Both distributions can be described by a Gaussian with the same 
fit parameters within errors (centroids 1.97±0.11 and 1.88±0.06 for 
the ULXs and the other sources, respectively, and widths 0.50±0.10 
and 0.41±0.05). The population of ULX candidates with power- 
law index T>3 is discussed in the text. 

This trend is illustrated in the upper panel of Figure 
where the ratio of the intrinsic luminosity derived using 
the mekal model to that derived using the power law 
model is shown against the power law index. Only a 
subset of ULX candidates are included in this figure: A 
sample with >500 detected counts and T<3 and all of the 
r>3 sources though, in both domains, only sources with 
statistically-acceptable fits from both the power law and 
the mekal model are included. The difference in derived 
intrinsic luminosities increases rapidly with increasing T. 
For the bulk of the sources modeled (those with r~1.8) 
however, both the mekal and power law models result in 
very similar intrinsic luminosities. 
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3 4 
Power Law Index 

Fig. 6. — (Upper panel) Ratio of the intrinsic luminosity derived 
using the mekal model to that derived using the power law model is 
shown against the power law index. The power law model predicts 
a higher intrinsic luminosity compared to the mekal model with the 
difference systematically increasing with power law index. (Lower 
panel) Corresponding mekal model plasma temperature decreases 
with increasing T. 
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This systematic effect helps explain the high-r tail ap- 
parent in the ULX candidate population of Figure but 
lacking in the lower luminosity source population. Ap- 
plying the mekal model to the 20 ULX and 7 non-ULX 
sources with T> 3 and acceptable power law fits resulted 
in mekal-derived intrinsic luminosities <10 39 ergs s _1 for 
17 of the 27 sources. These 17 are denoted in Table 2. 
Thus, adopting the mekal luminosities would result in 
similar high-r distributions for both populations. The 
lack, at the outset, of steep power law sources in the 
non-ULX population is due to the relatively low number 
of counts in the 0.5-8.0 keV band for this spectral shape 
combined with the fact that low-count sources were not 
selected for spectral analysis. However, the X-ray color 
analysis of § l3.2.2l shows that there are sources with colors 
equivalent to a steep power law in the low-count popula- 
tion. 

It is evident that more than one spectral model of- 
ten provides statistically-acceptable fits to the observed 
X-ray spectra. Often these models ostensibly represent 
quite different physical emission mechanisms. For exam- 
ple, physical models of X-ray phenomena that predict 
"soft" X-radiation, corresponding to a steep power law 
index, are cool thermal emission sources (lower panel, 
Figure 0) , such as supernovae, while phenomena that 
can be characterized as non-thermal, such as synchrotron 
emission from a power-law distribution of electrons or 
Comptonization, tend to result in relatively flat power 
law indices. Thus, while either model is statistically- 
acceptable, it is reasonable to assume that the majority 
of the steep power law sources are thermal sources. Many 
of these may be supernovae though three T>3 sources 
are among the most variable X-ray sources in the ULX 
sample (Table 2 and § I3.2.3JI . One of these is a known 
super-soft source in NGC 5457 (Mukai et al. 2003). 

Having established that both the power law and the 
mekal models predict similar luminosities for the major- 
ity of sources (those near r~1.8) and that both models 
effectively disriminate between "soft" and "hard" X-ray 
spectra, it remains to chose one or both models (or some 
other) as a baseline. The power law is arguably the eas- 
ier to grasp conceptually and requires the fewest number 
of parameters. Therefore, the absorbed power law model 
fit parameters and derived luminosities are reported in 
Table 2 whenever this model provided a statistically- 
acceptable fit. Additional models applied to sources not 
well-fit by an absorbed power law were thermal emission 
line models (mekal or vmekal), disk blackbody mod- 
els (diskbb) and combinations of these three (includ- 
ing power law) basic spectral models. Of the 24 ULX 
candidates poorly fit by the absorbed power-law model, 
11 were acceptably- fit using 1- or 2-component thermal 
emission line models, 9 using disk blackbody models, and 
4 were not well-fit by any of the trial spectral models. Pa- 
rameters for these models are reported in Table 2 (the 
power law parameters are quoted for the 4 spectra with 
no acceptable fits). 

Figure [3 shows the value of the power law spectral in- 
dex, r, against the number of detected source counts and 
Figure |H1 displays T against the derived intrinsic source 
luminosity. Not surprisingly, sources for which the simple 
absorbed power law model was unacceptable tend to be 
those sources with a high number of detected counts (e.g., 
a power law is unacceptable for 39% of the 51 sources 



with >1000 detected counts compared to only 8% of 177 
sources with <200 detected counts). 




100 1000 

Detected X— Ray Counts 

Fig. 7. — Results of simple absorbed power law model fits. 
Sources for which this model provides a statistically-acceptable 
fit at the 90% confidence level are marked with an X and those 
unacceptably-fit by a simple power law are denoted by triangles. 
Errors shown are extremes on the single interesting parameter 90% 
confidence intervals. They indicate typical values for all sources 
with power law models applied but are not displayed for sources 
with acceptable fits for clarity. Formal errors for the ULX candi- 
dates are provided in Table 2. Horizontal lines denote the average 
power law index, for sources with an acceptable fit, on domains 
spanning an equal number of sources. For clarity, these domains 
are separated by dotted vertical lines. 
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Fig. 8. — Results of simple absorbed power law model fits. Sym- 
bols are the same as in Figure |7| Note sources poorly-fit with 
a power law are uncorrelated with intrinsic luminosity and that 
many of the highest luminosity sources are well-fit with this simple 
model. Horizontal lines denote the average power law index, for 
sources with an acceptable fit, for ULX and for weaker sources. 
The vertical line denotes the estimated completeness limit of the 
sample of sources for which spectral modeling has been performed. 



The average value of the power law index, (r), for 
sources with statistically-acceptable power law fits, are 
1.79±0.09 (for sources with <137 c), 1.78±0.07 (137 to 
268 c), and 1.75±0.03 (>268 c) where the 3 ranges were 
chosen to contain equal numbers of sources. Similarly, 
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(r) is displayed in Figure for the ULX candidates 
(1.74±0.03), and for the weaker sources (1.77±0.04 above 
the completeness limit and 1.76±0.10 for sources below 
4 x 10 38 ergs s _1 ). The dispersion of T about the mean 
value is (7=1.18 for the ULX candidates and 0.43 for the 
non-ULX sources. The steep power law sources account 
for the larger dispersion in the ULX candidates popu- 
lation. Restricting the range to T<3 results in ct=0.51 
and 0.40 for the two groups. The average power law in- 
dex is therefore independent of the number of detected 
counts and of the intrinsic luminosities of the source 
above our completeness limit of ^4x 10 38 ergs s — 1 . Chan- 
dra studies of the nearby galaxies, e.g., M31 (Kong et al. 
2002) and M81 (Swartz et al. 2003), conclude that this 
trend continues ((F) ~ 1.8) even for sources as weak as 
3.4 x 10 37 ergs s _1 and 2 x 10 37 ergs s _1 , respectively. 

3.2.2. X-ray Colors 

There are many more weak sources detected than 
can be meaningfully analyzed through spectral fitting. 
A common tool used for low-count sources and low- 
resolution spectra is the X-ray color-color diagram. 
Prestwich et al. (2003) advocate {H - M)/(H + M + S) 
and (M — S)/(H + M + S) as colors that provide some 
physical insight into the nature of X-ray source popula- 
tions and these colors are adopted here. 

The X-ray colors of the 154 ULX candidates and 
of the 837 sources with estimated luminosities above 
10 38 ergs s^ 1 and >50 c are shown in Figures l9l and ITUl 
respectively. Error bars, propagated from the statisti- 
cal uncertainties in the three X-ray bands, are omitted 
from Figure ITUl for clarity. Colors of absorbed power law 
model spectra (for sources imaged on back-illuminated 
CCDs) are also shown for reference as solid curves. The 
front-illuminated devices have relatively lower response 
in the soft band so that similar curves for these devices 
do not extend to such low values of (M — S) /T. Never- 
theless, the differences are typically much less than the 
uncertainties in the X-ray colors. 
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Fig. 9. — X-ray colors of the ULX candidates. Solid curves 
denote colors of absorbed power law models of spectral indices 
r = 1, 2, 3, and 4 (from right to left) and for the range of absorbing 
columns nu = 10 20 to 10 24 cm -2 . Dashed curves denote constant 
absorption columns of nu = 10 20 , 10 21 , 2x 10 21 , and 5x 10 21 cm -2 
(from bottom to top). Errors shown were propagated from the 
statistical uncertainties in the three X-ray bands. 



- 


\ ''\ 

V \ 




1 






\ A/- y 














1 


\V • •" ■ 
\\ 






-0.5 


0.5 


1 






(H-M)/(H+M+S) 





Fig. 10. — X-ray colors of sources with luminosities, estimated 
from observed counts, in excess of 10 38 ergs s — 1 and detected 
counts >50. Curves same as Figure [5] Errors omitted for clar- 
ity. 



Highly-absorbed sources lie near (1,0) and super-soft 
sources near (0,-1) in these figures. The colors of steep 
power law sources are constrained along a narrow band 
extending from the super-soft source location at (0,-1) 
upwards and to the left towards ~ (-0.6,0.5). The region 
to the left of this band is physically inaccessible (curves 
of constant T converge rapidly as T increases as shown 
in the figures). According to Prestwich et al. (2003), 
the sources in the region centered at about (-0.2,0.0) are 
predominately low-mass X-ray binaries and are common 
to both spiral and elliptical galaxies. 

Visually, the two color distributions appear similar. 
Applying a two-dimensional Kolomogorov-Smirnov test 
(Peacock 1983; Fasano & Franceschini 1987) to the color- 
color distributions of the ULX candidates and of all 1900 
sources with luminosities above 10 38 ergs s _1 resulted 
in a significance level of Pks — 0.01. However, the 
Kolomogorov-Smirnov test does not account for (statis- 
tical) uncertainties in the data. When the test was ap- 
plied to the restricted set of sources with >50 detected 
counts (those depicted in Figures l9l and ITfj|> . the signif- 
icance level increased to Pks = 0.17, showing that the 
distributions are marginally consistent with being from 
the same parent population. 



3.2.3. X-ray Timing 

The Kolmogorov-Smirnov probability, Pks, that a 
source is constant was computed for the popula- 
tion of sources with estimated luminosities exceeding 
10 38 ergs s _1 and for the subset of ULX candidates. 
Analysis was limited to those sources with >50 detected 
counts as variability is more easily detected at higher 
signal-to-noise. 
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Fig. 11. — Results of a Kolmogorov-Smirnov test of the hypothe- 
sis that the distribution of photon arrival times for X-ray sources is 
equivalent to that of a source with constant flux. Shown is the value 
of the significance level, PfcSt f° r each source with estimated lumi- 
nosity exceeding 10 38 ergs s — 1 and >50 detected counts (dotted) 
and for the ULX candidates (solid) ordered by increasing values 
°f Pks- The line has slope unity. The abscissa is scaled by the 
inverse of the number of sources in the two populations (1/837 and 
1/154, respectively). 

Figure ITT1 displays the distribution of these probabili- 
ties for the two populations. The line denotes the linear 
trend of the data. For purposes of discussion, we may 
consider sources with Pffs<0.001 as significantly vari- 
able and those with Pffs<0.04, the point where the slope 
of the curve changes noticably, as likely to be variable. 
There are 80 sources (9.6%) with P KS < 0.04 and 33 
(3.9%) with P KS < 0.001 for the larger sample. If the 
ULX candidate population had the same distribution, 
then 15 and 6 ULXs, respectively, would be expected to 
lie below these values of Pks- For the ULX candidates, 
there are 22 and 7 sources, or 14.3% and 4.5%, respec- 
tively, below these two benchmarks. This is within sta- 
tistical errors of the expected values. Of particular im- 
portance is the fact that the photon arrival times for 86% 
of the ULX candidates are consistent with the constant 
flux hypothesis. 

Note that the ULXs tend to be those with higher num- 
bers of counts so that variability is more easily discerned 
with this test in this population. The K-S test is sensi- 
tive to long term variability even when the data has been 
finely binned. This is in contrast to other methods which 
often have an optimal binning typically of the same order 
as the variability timescale under consideration. Simu- 
lations were run to estimate this sensitivity. Consider 
the function R(l + Msin(27rf/P)), where R is the aver- 
age count rate, M is the modulation amplitude, and P 
is the period. Events were randomly generated in each 
timcbin with a Poisson distribution with the mean given 
by the function. For 5000 bins and R = 0.02 c bin -1 , 
the simulation will generate about 100 c. For this case 
and a period of 5000 bins, M needs to be about 0.35 in 
order for the K-S probability to be 10%. In other words, 
if the modulation is 35% and the K-S threshold is 10% 
then about half the time the K-S probability will be less 
than 10% and the variability "detected" . If P is 2500, 



then the amplitude needs to be 0.60 and for P of 1500 
an amplitude of 1.00 only results in a K-S probability 
of 15%. As expected the amplitude scales roughly as 
N~ x / 2 , where N is the number of counts. 

As displayed in Figure ITT1 the distributions are them- 
selves cumulative distribution functions and applying the 
Kolmogorov-Smirnov test to compare these two distribu- 
tions shows that they are consistent with being drawn 
from the same distribution (Pks = 0.95). 

3.2.4. X-ray Source Locations 

The relative radial position of each X-ray source can 
be expressed as the fraction, /, of the deprojected galaxy 
radius in units of 1/2 the D25 diameter. FigurelT2lshows 
the resulting surface distribution of all 3413 discrete X- 
ray sources detected in the survey, of all ULX candidates, 
of ULXs in spiral galaxies, and of ULXs in ellipticals. 
Recall that sources within 5"of the galaxy centers have 
been excluded from consideration as ULXs. 

The distributions can be approximated by a general- 
ized exponential function (plus a constant) of the form 

Aexp - ^/' 1 ) 11 " where h is a scale height and n is the 
index. The distributions of all detected sources and of 
all ULX candidates can be approximated by the same 
values of scale height (h = 0.06 ± 0.03) and index 
(n = 1.59 ± 0.32, where the quoted errors are for the 
larger population of all detected sources, x 2 = 105.9 for 
93 dof and 25.1 for 17 dof, respectively). Fits to the 
smaller subsets of spirals or ellipticals give large ranges 
to the fit parameter values when all parameters are al- 
lowed to vary. If n is held to the best-fit value of 1.59, 
then the values of the scale heights are within the uncer- 
tainties, h = 0.06 ± 0.04 (x 2 = 19.9 for 16 dof) for the 
spirals and h = 0.04 ± 0.02 (x 2 = 21.4 for 15 dof) for 
ellipticals. 
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Fig. 12. — Surface distribution of all X-ray sources detected in 
the sample (top) and of the subset of ULX candidates (bottom, 
full curve). ULX candidates in spiral galaxies are marked with 
triangles and those in elliptical galaxies with x's. The abscissa is 
the deprojected radial position, /, expressed as a fraction of the 
host galaxy's angular radius (= O.5D25) and the ordinate is the 
number of sources per unit area, fdf, on the range / to f + df. The 
curve drawn through the upper dataset is the best-fit generalized 
exponential function described in the text. 

Fitting a constant to the surface density distributions 
of the ULX populations in ellipticals and spirals results in 
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an acceptable fit for ellipticals but not for spirals. For el- 
lipticals, the value of the constant is 29.1 ±9.1, \ 2 = 23.3 
for 17 dof, on the entire domain. On the restricted do- 
main, / > 0.1, where the distribution appears to flatten 
(Figure EJ, the value is 28.8 ±9.1, x 2 = 13.9 for 15 dof. 
For spirals, \ 2 = 63.7 for 18 dof and \ 2 = 47.4 for 16 dof 
in these two domains, respectively. Thus, the ULX can- 
didates have the same spatial distribution, statistically, 
as does the larger population of all detected sources but 
the spatial distribution of ULX candidates in ellipticals 
is much flatter than for spirals with the exception of a 
few sources near the centers. 



3.2.5. Estimated Contribution from Background Sources 

The detection limit for discrete sources is about 3 x 
10~ 15 ergs cm~ 2 s _1 , for the average exposure time of 
39.6 ks, and the total area of the sky covered by the 
survey is 0.9 deg -2 . Comparing this to the Chandra 
Deep Field (CDF) results of Brandt et al. (2001) and of 
Rosati et al. (2002), we expect about 1100 of the 3400 de- 
tected sources are background (consistent with the value 
of the constant term, 1144 ±215, in the surface distribu- 
tion model described in H3.2.4|l . This corresponds to ~12 
background sources per galaxy given the average galaxy 
area is 36.2 min 2 . 

The results for the ULX candidates, however, are based 
on luminosities and not fluxes. Thus, a straightforward 
comparison to the CDF results is not possible. Instead, 
in FieurelLH is shown the \og(N(> S)) - log(S) distribu- 
tions of ULX candidates for each galaxy, where S is the 
observed flux and N(> S) is the number of sources with 
flux > S and scaled by the area of the galaxy within 
the instrumental field of view in units of deg -2 . Also 
shown is a curve, N(> S) = 16(S , /10 _13 ) _1 - 7 , represent- 
ing the high-flux log(iV(> S)) - log(5) results of the 
ASCA Large Sky Survey (Ueda et al. 1998) and Deep 
Sky Survey (Ogasaka et al. 1998) and the low-flux re- 
sults from Brandt et al. (2001). Points near or below this 
curve are likely background objects and not true ULXs. 
The number of background sources can be estimated as 
follows. If the flux from the weakest ULX candidate in a 
galaxy (of area A deg 2 ) is denoted S, then we expect a 
contribution of at most 16(S'/10 _13 ) _17 J 4 to the number 
of background sources from that galaxy field. If the slope 
of the galaxy's log(JV(> S)) - log(S) distribution is flat- 
ter than S -1 ' 7 , then this is an overestimate. The sum of 
such contributions from all galaxies gives a conservative 
estimate of 39 background sources among the 154 ULX 
candidates, or 25%. 

Note that the distribution depicted in Figure El sug- 
gests many ULX candidates in elliptical galaxies are 
more likely to be background sources than those in spiral 
galaxies. Using the method outlined above yields an esti- 
mate of 25 background sources, or 44%, among the ellip- 
tical galaxy ULX candidates while only 14 sources (14%) 
of the ULX candidates in spirals are potentially back- 
ground objects. This result is consistent with the sam- 
ple of elliptical galaxies analyzed by Irwin et al. (2003). 
However, there are several ULXs in our sample of ellip- 
tical galaxies that do not follow this trend. 
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Fig. 13.— log(A r (> 5)) - log(S) distribution for ULX candidates. 
The ordinate is the number of ULX candidates identified in each 
galaxy above a given flux scaled by the angular size of that galaxy. 
Spiral galaxies are denoted by triangles and solid lines. Ellipticals 
are denoted by crosses and dotted lines. The power law curve 
denotes an approximate fit to the point-source X-ray background 
(e.g., Brandt et al. 2001), scaled to our bandpass (0.5-8.0 keV), 
and is given by JV(> S) = 16(S/10- 13 ) -1 - 7 . Note several galaxies 
host only one ULX candidate. 

3.2.6. X-ray Luminosity Function 

The cumulative luminosity functions (XLFs) of the 
ULX candidates in elliptical and in spiral host galaxies 
are displayed in Figure El There are 57 ULX candi- 
dates in the 27 elliptical galaxies and 97 in the 55 spi- 
ral galaxies. While the number per galaxy are roughly 
equal, ULXs in spiral galaxies are much more luminous. 
Two-thirds of the ULXs in spiral galaxies are above 
2 x 10 39 ergs s -1 while only 1/3 of ULXs in ellipticals 
are this luminous. 
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Fig. 14. — Cumulative luminosity function of ULX candidates 
in spiral (upper curve) and in elliptical galaxies. A power law with 
slope -1.72±0.14 describes the luminosity function for ellipticals. 
A broken power law model is shown for spirals. (A power law with 
an exponential cut-off also provides an acceptable fit for spirals 
though a single power law does not.) The slope is -0.56±0.05 at 
low luminosity and the break luminosity is ~10 X 10 39 ergs s _1 . 

Using a maximum-likelihood statistic, a power law ad- 
equately describes the luminosity function for ellipti- 
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cals with a slope -1.72±0.14. More complex models do 
not improve the fit significantly. A broken power law 
(change in statistic AC = 26.7 compared to a simple 
power law) or a power law with an exponential cut-off 
(AC = 23.8) is needed to describe the luminosity func- 
tion for spirals. The slope of the broken power law 
is -0.56±0.05 at low luminosity, the break luminosity 
is 9.6±0.6 in units of 10 39 ergs s^ 1 , and the slope at 
high luminosity is -1.88±0.30. The slope of the cut-off 
power law model is -0.64±0.09 with a cut-off luminosity 
of (19.4 ± 8.6) x 10 39 ergs s" 1 . 



3.3. ULXs and their Host Galaxies 

Several global prop ertie s of the host galaxies are avail- 
able as described in § 12.21 Linear correlation coefficients 
were computed for various paired combinations of the fol- 
lowing quantitites compiled for each galaxy: number of 
ULXs (iVuLx); total X-ray luminosity from ULXs (Ltjlx) 
average ULX luminosity (£tjlx/Atjlx); blue luminosity 
{Lb), far-infrared luminosity (£fir), nearest-neighbor 
distance (6 P ), and D25 area. Linear correlation coeffi- 
cients for the most significant correlations are listed in 
Table 3. In addition to the entire sample of galaxies, 
coefficients were computed for the subsets of elliptical 
and of spiral galaxies and for the subsample of Chandra 
galaxies with the same distribution of absolute B mag- 
nitudes as the Tully catalogue of galaxies (see § 13.1(1 . 

The strongest correlations for elliptical galaxies are 
between ULX properties and the host galaxy's Lb, 
while those for spirals are with Lfir and a marginally- 
significant correlation with 6 p . There are no significant 
correlations for the subsample of elliptical galaxies. The 
subsample of spiral galaxies again correlate with Lfir 
but also with Lb- The subsample of galaxies was de- 
signed to minimize some of the biases in the full sample. 
Although the subsample contains only 34 galaxies, the 
significance of the correlations betwen ULX properties 
and L-piB, for spirals remains strong enough to show this 
result is robust. 

While there are ^2 ULX candidates per galaxy for 
both ellipticals and spirals 1 ^3.2.6(1 . the number per unit 
(10 42 ergs s _1 ) blue luminosity for ellipticals is only 
0.11±0.02 compared to 0.30±0.11 for spirals. The total 
ULX luminosity, Ltjlx, per unit Lb is 0.19±0.05 for el- 
lipticals and 1.63±0.88 for spirals. Therefore, the pairs of 
parameters were again examined after first normalizing 
to unit Lb- No significant changes in the correlation co- 
efficients occurred though the correlations between ULX 
properties and Lfir increased for spiral galaxies. 

The correlation between the number of ULX candi- 
dates and Lfir is displayed in Figure for the spi- 
rals and ellipticals. For spirals, the linear relation be- 
tween the two is iV UL x=(0.022±0.01)L F iR/10 42 -l-(0.64± 
0.3). Note that the highest Lfir bin is dominated 
by the Antennae pair of galaxies with 13 ULX can- 
didates. On the other hand, the spiral galaxies with 
£fir^4 x 10 42 ergs s _1 are the dwarf irregulars (cf. Fig- 
ure^). For ellipticals, the strongest correlation is be- 
tween the number of ULXs and blue luminosity and is 
given by ^ UL x=(0.095 ± 0.05)L B /10 42 + (-0.11 ± 0.7). 
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Fig. 15. — Number of ULX candidates per galaxy is shown 
against the host galaxy's far-infrared luminosity. For spiral galaxies 
(solid histogram), the correlation is significant (Pearson's linear 
correlation coefficient r = 0.88) but not for ellipticals (dotted, r = 
0.15). The numbers of elliptical/spiral galaxies within each LpiR 
bin are displayed across the top of the figure. 

4. DISCUSSION 

We have presented the spectrophotometric X-ray prop- 
erties of a population of 154 ULX candidates taken from 
a sample of 82 galaxies observed with Chandra. We 
have shown that ULXs are more numerous and luminous 
in galaxies with indications of recent star formation as 
measured by their far-infrared luminosities and by their 
morphologies (Hubble type and, less significantly, near- 
est neighbor distance) . We have also compared the ULX 
candidates with the less-luminous X-ray sources in the 
same galaxy fields and found their X-ray properties sta- 
tistically indistinguishable. 

4.1. ULXs and Less-Luminous Sources 

The similarity between ULX candidates and the less- 
luminous population exists among all observable X-ray 
properties: spectral shapes, colors, time series, and (ra- 
dial) positions within their host galaxies. This result has 
two important implications. First, as the less-luminous 
population is a composite of low- and high-mass X-ray 
binaries and supernova remnants (to our formal detec- 
tion limit but also must include some foreground objects 
such as cataclysmic variables and stellar coronae as well 
as background AGNs); our results imply the ULX can- 
didates are also a heterogeneous group of objects. Sec- 
ondly, the similar distributions of physical characteristics 
suggests, as the simplest explanation, that ULXs may be 
the high-luminosity end of a continuous distribution of 
supernova remnants and accreting systems such as X-ray 
binaries. 

These conclusions do not diminish the fact that ULXs 
retain their distinction as the most luminous non-nuclear 
X-ray sources. How are such high luminosities achieved? 
Contemporary explanations for the ULX phenonema fall 
into three broad catagories; massive accreting objects, 
anisotropically-emitting sources, and supernovae. As 
part of their appeal, each of these models effectively 
circumvent the benchmark Eddington limit luminosity, 
L~1.4x 10 38 (M/M Q ) ergs s _1 , that expresses the balance 
between the inward force of accreting material against 
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the outward radiation force. X-ray binaries with com- 
pact object masses >100 Mq raise the Eddington limit 
to >10 40 ergs s _1 . Anisotropically-emitting X-ray bina- 
ries avoid the Eddington limit in two ways: The emission 
may be in a direction other than opposing accretion and, 
secondly, the actual luminosity is less than the appar- 
ent (assumed isotropic) luminosity by a factor b = £l/Air 
where is the solid angle of emission. For accreting 
sources, it is also necessary to attain an accretion rate, 
A/>10 _6 6 M yr _1 , as the apparent luminosity scales as 
~6x 10 45 M/6 ergs s _1 [and so the slim-disk scenario that 
applies at high accretion rates (Abramowicz et al. 1988; 
Watarai et al. 2001; Ebisawa et al. 2003) may be favored 
over standard thin-disk models]. The third catagory, su- 
pernovae, are, of course, not Eddington limited in the 
first place. Each of these scenarios predict different ob- 
servable X-ray properties which can be compared to our 
results. 

4.1.1. Intermediate- Mass Black Holes 

Colbert & Mushotzky (1999; see also Makishima et 
al. 2000; Colbert & Ptak 2002) suggest that ULXs are 
accreting black holes (BHs) with masses intermediate be- 
tween stellar mass (^20 M©) and supermassive objects 
(active galactic nuclei, 10 6 to 10 9 M Q ). This is perhaps 
the most intriguing scenario because it suggests a possi- 
ble observable link between stellar collapse and the for- 
mation of active galactic nuclei (see van der Marel 2004; 
Miller & Colbert 2003 for reviews). Fryer & Kalogera 
(2001) show, however, that single stars in the current 
epoch rarely form BHs with masses >15 Mq. First- 
generation, zero-metallicity (Population III) stars, on 
the other hand, were on average more massive than 
today's population, would have retained more mass as 
they evolved because of inefficient wind mass loss, and 
may have left more massive BH remnants. Intermedi- 
ate mass black holes (IMBHs) may also form in dense 
young star clusters through runaway merging of massive 
stars (Ebisuzaki et al. 2001; Portegies Zwart & McMillian 
2002; Portegies Zwart et al. 2004) or through the grad- 
ual accural of compact stellar remnants onto a seed black 
hole over the lifetime of the cluster (Miller & Hamilton 
2002). 

In analogy to stellar-mass black hole X-ray binaries 
(e.g., McClintock & Remillard 2003), the IMBH X-ray 
spectrum should be a power law with a soft thermal 
(blackbody) component with a characteristic tempera- 
ture scaling as the -l/4th power of the compact object 
mass (e.g., Makishima et al. 2000). Those ULX candi- 
dates with blackbody disk spectral shapes typically have 
temperatures of order 1 keV (with two exceptions which 
lie near (0,-1) in the color-color diagram, Figure |§1 i.e., 
they are supersoft sources with no additional power law 
component). Thus, we find no prominent component at 
fcT~0.1-0.2 keV in our sample. A few ULX candidates 
with soft accretion disk components have been reported 
(e.g., Miller et al. 2003; Roberts & Colbert 2003). 

To further the analogy to stellar-mass BHs, IMBHs 
can also be expected to display soft/hard spectral tran- 
sitions. Furthermore, transient behavior is predicted by 
Kalogera et al. (2004) for accreting IMBH systems due 
to a thermal-viscous disk instability. Neither of these 
transitions are accessible in the short observations typi- 



cal of our sample but have been reported for a few well- 
monitored ULXs (e.g., La Parola et al. 2001; Kubota et 
al. 2001). 

IMBHs formed from individual Population III stars are 
expected to be distributed throughout the host galaxy 
and to be present in all galaxy types (Madau & Rees 
2001). IMBHs formed in dense stellar environments 
should remain associated with globular clusters. They 
should occur frequently in elliptical galaxies as these 
galaxies have large numbers of globular clusters. Further- 
more, the ULX spatial distribution should follow that of 
the halo stars in contrast to the distribution of weaker 
sources that follow the distribution of disk and bulge 
stars. These trends are difficult to discern in the current 
sample. The ULXs are distributed in the same way as 
the weaker sources (though the distribution in ellipticals 
is arguably flatter, ^3.2.41) and the number of ULXs per 
galaxy is no higher in cluster-rich ellipticals compared 
to spirals. Note that an IMBH accretor need not be 
invoked for the majority of ULXs in elliptical galaxies 
where 2/3 of the ULXs f§ I3.2.6[) are less luminous than 
2 x 10 39 ergs s _1 and, according to Irwin, Bregman, & 
Athey (2004), those more luminous than 2x 10 39 ergs s _1 
may all be background objects. 

4.1.2. Anisotropic X-ray Binaries 

There are several beaming models: One model (King 
et al. 2001) assumes a phase of super-Eddington accre- 
tion is present, perhaps due to thermal-timescale mass 
transfer, that forms a thick disk with a central funnel 
that results in mildly-beamed emission (see also Madau 
1988). Another model envisions the formation of a jet 
with an enhancement of the X-radiation via Compton 
scattering of photons (from the accretion disk or from a 
high-mass companion) by relativistic electrons in the jet 
(Georganopoulos et al. 2002). A third model envisions 
direct, relativistic Doppler boosted, synchrotron emis- 
sion from the (face-on) jet itself (Fabrika & Mescheryakov 
2001). 

The spectral signature of the beaming model has not 
been specified. Compton scattering from jets and di- 
rect synchrotron emission should display a characteristic 
steep power law spectra. If ULX emission is beamed 
and weaker sources are not, then one would expect the 
two populations to have different spectral indices. While 
the spectral shape of the majority of the ULX candi- 
dates are consistent with power laws, their spectral in- 
dices do not differ from the less luminous population. 
Large-amplitude variability is also expected from rela- 
tivistic jets. The high mass transfer rates implied for 
the mildly beamed sources (King et al. 2001; Kalogera 
et al. 2004) result in stable disks and persistent X-ray 
emission. There are no preferred locations or galaxy 
type for beaming models. However, thermal-timescale 
mass transfer occurs when the main-sequence mass-donor 
(companion) star is more massive than the compact ac- 
cretor. This requires an early-type companion of mod- 
erate mass, ^5 M Q , and hence such systems will occur 
rarely in ellipticals. 

4.1.3. Energetic Supernovae 

A third scenario is an unusually energetic supernova, or 
hypernova (Wang 1999). Supernovae exploding in high- 
density environments have been hypothesized to reach 
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X-ray luminosities sufficient to power the broad-line re- 
gions of active galactic nuclei (Terlevich, Melnick, & 
Moles 1987; Terlevich et a!. 1992). Supernovae have been 
observed with luminosities as high as ^10 40 ergs s -1 
months to years after their discovery (Schlegel 1995; 
Immler & Lewin 2003). For example, SN 1986 J in 
NGC 891 has a luminosity of ~6 x 10 39 ergs s _1 accoding 
to our calculations (Table 2). 

Supernovae display the distinct thermal emission-line 
X-ray spectrum from shock-heated plasma; often with a 
harder bremsstrahlung component from a forward shock. 
This spectrum is easily distinguished from a power law 
shape. Only ~7% of ULX candidates require a t hermal 
model to reproduce their spectra. As shown in § 13.2.11 
about 15% of the ULX candidates with acceptable power 
law fits have steep slopes, r>3, suggesting a thermal 
spectrum. Thus, at least 22% of the ULX candidates 
are likely thermal sources. The actual number could be 
higher since a thermal model does provide a formally- 
acceptable fit to 78% of the ULX candidates. Neverthe- 
less, it is unlikely that the majority of ULX candidates 
are thermal sources since the average power law model 
slope of 1.8 (or the average thermal model temperature 
of 18.4 keV) is much more descriptive of non-thermal 
sources. 

Supernovae should be steady sources on short 
timescales and decline slowly over periods of months to 
years. Again, while results of variability tests are consis- 
tent with the majority of ULX candidates being steady 
sources, individual observations are too short to discern 
long-term trends and some low-level variability may be 
hidden in the low-count light curves. Variability has been 
reported for several ULX candidates based on multi- 
epoch monitoring (e.g., Fabbiano et al. 2003; Roberts et 
al. 2004, La Parola et al. 2001). In combination with de- 
tailed spectral modeling and radio wavelength imaging, 
the number of energetic supernovae contributing to the 
ULX population should soon be much better constrained. 
Core collapse supernovae have massive star progenitors 
consistent with the association of ULXs with late-type 
star-forming galaxies. 

4.2. ULXs and their Host Galaxies 

It has long been known (David, Jones, & Forman 1992) 
and recently confirmed (Ranalli, Comastri, & Setti 2002; 
Grimm, Gilfanov, & Sunyaev 2003) that the total X-ray 
luminosity of a galaxy correlates with recent star forma- 
tion rate as measured by its far-infrared luminosity, Z/fir- 
We have shown that this trend extends to the numbers 
and total luminosities of ULXs found in such galaxies. 
This result was anticipated by the large numbers of ULXs 
in some FIR-bright galaxies (e.g., NGC 4038/4039, Zezas 
et al. 2002; NGC 3256 Lira et al. 2002; NGC 4485/4490, 
Roberts et al. 2002) and for a small sample of galaxies 
by Grimm, Gilfanov, & Sunyaev (2003). A similar con- 
clusion can be indirectly deduced through studies of the 
cumulative X-ray luminosity functions of nearby galax- 
ies: The flat slopes of the XLFs implies that the most 
luminous sources dominate the total X-ray luminosity of 
spiral and starburst galaxies (e.g., Kilgard et al. 2002; 
Grimm et al. 2003; Colbert et al. 2004). 

Two of the theoretical explanations for ULXs discussed 
in the previous section predict an association of ULXs 
with recent star formation. These are the mildly-beamed 



thermal-timescale mass transfer X-ray binaries described 
by King et al. (2001) and X-ray bright supernovae. While 
a case might be envisioned for IMBHs formed in prompt 
collapse of the cores of young super-star clusters, studies 
of the Antennae by Zezas et al. (2002, see also Clark et al. 
2003) and of three other starburst galaxies by Kaaret et 
al. (2003) determined ULXs are often near but not spa- 
tially coincident with these clusters as would be expected 
if they are IMBH binaries. Recent numerical simulations 
(Portegies Zwart et al. 2004) indicate the formation fre- 
quency of IMBHs in super-star clusters depends sensi- 
tively on cluster initial conditions and dynamical fric- 
tion timescale and thus does not simply scale with the 
host galaxy's (global) star formation rate. As only a 
small fraction of ULX candidates in our survey display 
the requisite thermal emission-line spectrum characteris- 
tic of supernovae, there remains only the model of King 
et al. (2001) as a reasonable representation of the major- 
ity of ULXs. 

Does the model of King et al. (2001) provide a viable 
explanation for the results of the present survey? The 
model is essentially that of a high-mass X-ray binary 
(HMXRB) undergoing an episode of thermal-timescale 
mass transfer. Such systems arise when the donor star 
envelope is in a radiative phase either when it is more 
massive than its companion or first fills its Roche lobe 
while expanding to the red giant stage. Particulars of 
the mass transfer mechanism, possibilities of common 
envelope formation, and details of beaming and resulting 
luminosities are beyond the scope of the present work. 
However, the question remains if enough systems of this 
type occur to account for the numbers of ULXs reported 
here. 

Assuming ~2% of OB stars form HMXRBs (e.g. Dal- 
ton & Sarazin 1995, Helfand & Moran 2001), that the 
thermal-timescale phase always occurs and lasts ^10 5 yr 
(King et al. 2001) or 1% of the lifetime of a star (10-40 
Myr for stars of initial mass 8-20 M Q ), and a beaming 
factor b = 0.1; then at least ~5 x 10 4 OB stars per ULX 
are required. From results of evolutionary synthesis mod- 
els, Leitherer & Heckman (1995) find this many O stars 
are formed for a star formation rate (SFR) of ~2 M© 
yr -1 over a period of >10 Myr. This result is for a par- 
ticular choice of initial mass function, metallicity, upper 
mass cutoff, and star formation time-scale representative 
of conditions found in infrared-luminous starburst galax- 
ies. The relation between Lfir and SFR deduced by Lei- 
therer & Heckman (1995) is 9 SFR=0.045(L F ir/10 42 ) M 
yr -1 (see also Kennicutt 1998). These values predict ^2 
to 5 ULXs per Lfir=10 44 ergs s _1 , though the uncer- 
tainties are large. From H3.3I the number of ULXs per 
(spiral) galaxy is ~1 for galaxies with Lfir~10 43 ergs s _1 
rising to ^5 for Lfir^IO 44 ergs s _1 . Thus, the number of 
ULXs in high Lfir galaxies follows the trend expected for 
a short-lived high stellar mass origin. Leitherer & Heck- 
man (1995) also predict a peak massive-star supernova 
rate of 0.02 yr^ 1 at a SFR of 1 M Q yr" 1 which will also 
contribute to the ULX population of active star-forming 
galaxies in proportion to £fir- 

For lower £fir, and particularly for elliptical galax- 

9 I/FIR m this expression is the 8-1000 (im luminosity which is 
about 50-100% larger than the 42.5—122.5 fim luminosities listed 
in Table 2 assuming a warm dust temperature. 
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ies, the number of ULXs is larger than this simple esti- 
mate predicts. Being dominated by the elliptical galax- 
ies, ULXs in this group are correlated most strongly with 
the host galaxy's Lb- The number of ULX candidates per 
elliptical galaxy scales roughly as Aulx^0.1(Lb/10 42 ) 
(§ 13.3(1 . It is thus natural to relate the lower-luminosity 
ULXs to the long-lived population of stars, namely, to 
the low-mass X-ray binaries (LMXRBs). ULX candi- 
dates in weak Lfir galaxies are also less-luminous than 
those in active star-forming galaxies as the average ULX 
luminosity also correlates with Lfir (§ 13.3(1 . The XLFs 
depicted in Figure ITU show that the ULXs in ellipticals 
are less-luminous than those in spirals. Irwin, Athey, & 
Bregman (2003) and Irwin, Bregman, & Athey (2004) 
show further that all ULX candidates in elliptical galax- 
ies more luminous than 2 x 10 39 ergs s^ 1 are likely back- 
ground sources. 

LMXRBs with stellar-mass (Af-10-15 M Q ) BH ac- 
cretors can reach luminosities of ~2 x 10 39 ergs s _1 
without violating the Eddington limit. The difficulty 
is achieving accretion rates of up to ^10" 6 M© yr" 1 . 
Although such high rates are seen in the soft X-ray tran- 
sients (King 2002; Terashima & Wilson 2004), persistent 
sources within our Galaxy do not achieve this rate. A 
typical Galactic LMXRB radiating at a few 10 37 ergs s _1 
has an accretion rate of ^10~ 8 M yr" 1 and can main- 
tain this luminosity for M/M~10 8 yr for a 1 M Q mass- 
donor star. If all accretion rates are equally likely, then 
the probability of observing a ULX in a single observation 
would be inversely proportional to the lifetime at that 
accretion rate (e.g., Wu 2001). For this assumption, a 
luminosity of ~10 37 ergs s _1 is expected to be 100 times 
more common than ~10 39 ergs s _1 . This is consistent 
with what is inferred from the ULX cumulative luminos- 
ity function for ellipticals, Figure El There is one ULX 
per 100 X-ray sources above 6.3x 10 37 ergs s^ 1 . Of order 
100 X-ray sources per Lb=10 43 ergs s _1 are routinely 
detected above a few 10 37 ergs s _1 in elliptical galaxies 
(e.g., NGC 4697, Sarazin et al. 2001) and in the bulges 
of spiral galaxies (e.g., NGC 3031, Swartz et al. 2003). 
The average blue luminosity of the elliptical galaxies in 
our sample is (18 ± 3) x 10 42 (Figure EJ|. Thus, the ob- 
served rate of ~2 ULX candidates per elliptical galaxy 
(or 1 ULX candidate per Lb = 10 43 ergs s _1 ) is consistent 
with a LMXRB origin. The actual number is somewhat 
less since nearly 50% of the ULX candidates in elliptical 
galaxies are potentially background sources (§ 13.2.5(1 but 
this line of reasoning suggests that the dominant factor 
determining the slope of the cumulative X-ray luminos- 
ity function for ellipticals is simply the lifetimes of the 
LMXRBs. 

Low-mass XRBs also exist in spiral galaxies. However, 
inspection of Figure0|shows that, for spirals, Lb is more 
a measure of the young stellar content than of the total 
number of stars (or galaxy mass). Thus, while the con- 
tribution of LMXRBs to the ULX population is expected 
to scale with galaxy mass, it may not vary linearly with 
Lb for spirals in the present sample. Colbert et al. (2004) 
used B and K s luminosities as a measure of galaxy mass 
to conclude that more than 20% of ULXs in spiral galax- 
ies originate from the older stellar population. 

Arguments for contributions from the young and the 
old stellar populations to the total X-ray luminosity of 



nearby galaxies have been made since the Einstein Obser- 
vatory era (see the review by Fabbiano 1989). More re- 
cently, Colbert et al. (2004) show that X-ray point source 
populations can also be described as a superposition of 
contributions from these two stellar populations. Here it 
has been shown that the ULX candidate population also 
displays this dichotomy. Considering ULXs as a subset 
of all X-ray sources or, equivalently, as a subset of the 
total X-ray luminosity of a galaxy, this result supports 
the idea that ULXs are the high-luminosity end of a dis- 
tribution of HMXRBs and supernovae and of LMXRBs. 

4.3. Summary 

X-ray properties and celestial positions of a sample 
of 154 ULX candidates have been tabulated. The X- 
ray properties alone are unable to discriminate between 
ULXs and the less-lumious population (to a formal com- 
pleteness limit of ~4x 10 38 ergs s" 1 ). This could perhaps 
be foreseen since, for example, the power law indices of 
accreting X-ray sources ranging from stellar-mass neu- 
tron stars and black holes to supermassive active galac- 
tic nuclei are, on average, comparable to the mean power 
law index, r=1.74±0.03, obtained here for ULX candi- 
dates. Based on spectral shape, namely a steep power 
law index or, equivalently, a cool plasma temperature, 
about 20% of the ULX candidates are potentially super- 
novae but this group also includes a few variable and/or 
supersoft sources. The supersoft (or 'quasisoft') spectra 
may be a signature of IMBHs as these objects are ex- 
pected to have a soft blackbody disk component. Most 
of the sources studied here, however, can be described 
by a simple power law and do not require an additional 
(soft) thermal component. Future spectrophotometric 
X-ray observations will help distinguish between tran- 
sient XRBs and slowly-evolving young supernovae but 
optical and radio data (and rapidly-sampled X-ray tim- 
ing data) will be needed to discern, e.g., beamed sources 
from IMBHs. The precise source locations provided by 
Chandra should be a tremendous aid to such followup 
observations. 

The number of ULX candidates per galaxy is roughly 
the same for spirals and ellipticals but the luminosities 
of the ULXs are ~10 times higher in spirals compared to 
ellipticals and the ULX candidate population in ellipti- 
cals is severely compromised by background sources (esti- 
mated to be 44% in the present sample compared to only 
14% for spirals). There is also a strong correlation be- 
tween the number and average luminosity of ULXs in spi- 
rals and their host galaxy's far-infrared luminosity; and a 
weaker (anti-)correlation with nearest-neighbor distance. 

The X-ray data do not favor any particular model for 
ULX phenomena. The similarity to properties of weaker 
sources strongly suggests ULXs are not a distinct class 
but composed of a heterogeneous mixture as are the 
weaker sources; in which case all models remain viable 
on a case by case basis. The association with recent 
star formation, and the significant lack of bright ULXs 
in elliptical galaxies, does make a strong case for young, 
short-lived, systems such as HMXRBs as the dominant 
component of the more luminous ULX population. On 
the other hand, the population of ULXs in ellipticals and 
at least some of the ULXs in spirals can be explained as 
the high-luminosity end of the LMXRB population. 

The intrinsic luminosities of the ULX candidates in the 
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sample range up to ~10 x 10 39 ergs s _1 in the 0.5-8.0 keV 
band, the point where the luminosity function rolls over, 
with a few objects radiating at even higher luminosities. 
The Eddington luminosity limit thus requires accretor 
masses of order 80 M and possibly even higher. Stellar- 
mass BHs can form with masses up to ~20 M Q (Fryer & 
Kalogera 2001). If a companion of high initial mass can 
add some 10-20 M Q during a phase of super-Eddington 
accretion and the system is observed near the endpoint 
of this phase, then the luminosity of such a high mass 
X-ray binary could readily be of order 10 x 10 39 ergs s _1 
if the Eddington limit is exceeded by modest amounts 
and a high accretion rate can be acheived. While it may 
still be difficult to attain the highest luminosities inferred 
in the current study, this need only occur rarely to ac- 
count for the handful of very luminous ULXs (among the 
thousands of sources detected) . 
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0.67 


3.1 


41 
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0.40 


0.35 
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12.9 


42 
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24.32 
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0.95 


19.9 


43 
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6.82 
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2.2 
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2000-04-29 


51.0 


0.79 


1.7 


44 


NGC 5253 


Pec 


3.2 C 


0.51 


1.65 


10.7 


2032 


2001-01-13 


56.6 


1.00 


1.2 




NGC 5457 
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7.0 C 


19.98 


35.74 


10.1 


934 


2000-03-26 


98.2 


0.53 


2.9 


45 


Circinus 
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4.0 v 


4.94 


23.02 


3.1 


356 


2000-03-14 


23.1 


1.00 


9.8 


46 


NGC 5774 
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26.8 v 


4.86 


4.54 


1.0 


2940 


2002-04-05 


58.2 


1.00 


76.2 




NGC 5775 


SBc 


26.7 v 


17.96 


120.51 


1.0 


2940 


2002-04-05 


58.2 


1.00 


75.6 




NGC 6503 


SA(s)cd 


4.1 t 


1.18 


1.39 


41.0 


872 


2000-03-23 


13.2 


1.00 


8.0 




NGC 6946 


SAR(rs)cd 


5.5 t 


18.12 


31.53 


25.3 


1043 


2001-09-07 


58.3 


0.87 


4.8 


47 


NGC 7320 


SA(s)d 


13.8 v 


1.41 


1.37 


13.9 


789 


2000-07-09 


19.7 


1.00 


68.3 


48 


NGC 7331 


SA(s)b 


14.6 C 


29.33 


65.75 


27.7 


2198 


2001-01-27 


29.5 


0.62(1/2) 


51.4 




IC 1459 


E3-4 


29.2 I 


18.43 


1.85 


1.2 


2196 


2001-08-12 


58.8 


1.00(1/2) 


49.8 


49 


IC 5332 


SA(s)d 


9.4 v 


2.94 


0.96 


14.4 


2066 


2001-05-02 


52.1 


1.00 


9.8 


13 



References. — (1) Strickland et al. 2000; (2) Eraclcous ct al. 2002; (3) McDowell et al. 2000; (4) Jeltema et al. 2003; (5) Brcgman & Irwin 
2002; (6) Smith & Wilson 2003: (7) Irwin, Sarazin & Bregman 2002; (8) Kim & Fabbiano 2003; (9) Irwin, Athey & Bregman 2003; (10) Blanton, 
Sarazin & Irwin 2001; (11) Hcikc ct al. 2003; (12) Schlegel & Pannuti 2003; (13) Humphrey ct al. 2003; (14) Swartz ct al. 2003; (15) Matsumoto 
ct al. 2001; (16) Cecil, Bland-Hawthorn & Veilleux 2002; (17) Ott, Martin & Walter 2003; (18) Schlegel 2001; (19) Wang, Chaves & Irwin 2003; 
(20) Strickland et al. 2001; (21) Fabbiano, Zczas & Murray 2001; (22) Collingc et al. 2001; (23) Yang, Wilson & Ferruit 2001; (24) Cagnoni et al. 
2003; (25) Wilson, Yang & Cecil 2001; (26) Jimcz-Bailn ct al. 2003; (27) Kundu ct al. 2003; (28) Finoguenov & Jones 2002; (29) Sivakoff, Sarazin 
& Irwin 2003; (30) Iwasawa ct al. 2003; (31) Summers ct al. 2003; (32) Maccaronc, Kundu & Zepf 2003; (33) Roberts ct al. 2002; (34) Wilson 
& Yang 2002; (35) Di Stefano et al. 2003; (36) Wang ct al. 2001; (37) Jones et al. 2002; (38) Sarazin, Irwin & Bregman 2001; (39) Pellegrini 
ct al. 2002; (40) Schurch, Roberts & Warwick 2002; (41) Kraft et al. 2001; (42) Roberts et al. 2001; (43) Liu ct al. 2002; (44) Soria & Wu 2002; 
(45) Mukai ct al. 2003; (46) Smith & Wilson 2001; (47) Holt et al. 2003; (48) Trinchicri ct al. 2003; (49) Fabbiano ct al. 2003 

a C:Cepheid period-luminosity; I:I-band SBF; T:Tip of RGB; P:Pla notary nebula luminosity function; G:Globular cluster luminosity function; 
t:ter tiary estimators; v:Tully-Fisher and variants including Tully (1988); F:averagc of primary estimators from Fcrrerasc ct al. (2000) 
b Nearest-neighbor distance in units of the -D25 diameter 
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TABLE 2 
Properties of ULX Candidates 
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6 
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8 
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4 
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8 
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29 


5 
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8 
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3 
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.4 
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10 


53. 


2 
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9 
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4 
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3079 


10 


02 


02 
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38 
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.8 
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10 
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27 
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,7 
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10 
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34 


56. 


9 
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1 
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.8 
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16 
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03 
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1 
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7 
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11 
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12 


01 


51 
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01 


51 
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52 
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01 


52 
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6 
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01 


52 
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52 
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.8 


NGC 


4038 
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01 


53 
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Spectral™ 
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(10 39 erg s 1 ) 


(10 39 erg s _1 ) 


color 


color 
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1.50 + 0.08 


1.78 + 0.08 


-0.17 


-0.12 


0.001 


3.19/5 


0.99 + 0.22 


1.08 + 0.24 


-0.05 


-0.14 


0.521 


4.38/3 


1.16 + 0.24 


1.57 + 0.52 


-0.13 


0.28 
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1.07 + 0.19 


1.46 + 0.34 


-0.25 


-0.09 
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3.55/4 
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1.32 + 0.22 
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0.16 


0.274 
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4.05 + 0.26 
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-0.22 


0.40 
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79.88/77 


5.03 + 0.16 


8.33 + 0.35 


-0.08 


0.42 


0.603 


75.10/61 


11.33 + 0.42 


13.51 + 0.47 


0.28 


0.31 


0.822 


20.82/20 


1.43 + 0.09 


1.49 + 0.11 


-0.25 


0.09 


0.198 


12.20/16 


0.63 + 0.08 


3.07 + 1.35 


-0.39 


0.02 


0.613 


20.78/20 


0.62 ± 0.04 


1.21 + 0.16 


-0.30 


-0.09 


0.715 


19.20/20 


0.84 + 0.06 


1.17 + 0.10 


-0.28 


-0.01 


0.700 


5.35/9 


0.52 + 0.12 


1.11 + 0.54 


-0.32 


0.11 


0.766 


8.50/8 
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3.44 + 1.57 


-0.23 


-0.30 
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1.10/3 


1.30 + 0.24 


1.39 + 0.24 


-0.02 


-0.08 


0.371 


8.62/9 


1.11 + 0.15 


1.31 + 0.16 


-0.18 


-0.22 


0.351 


4.83/6 


1.36 + 0.19 


1.45 + 0.17 


-0.13 


0.02 


0.268 


17.59/15 


1.07 + 0.10 


2.26 + 0.58 


-0.26 


-0.37 


0.022 


23.76/25 


3.53 + 0.22 


3.67 + 0.24 


-0.15 


-0.07 


0.622 


11.84/7 
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0.83 + 0.41 


-0.39 


0.24 
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5.84/8 
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-0.28 


-0.28 
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15.93/12 
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-0.13 


0.306 


9.06/6 


1.17 + 0.21 


1.74 + 0.42 


-0.27 


0.30 
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4.82/8 


1.60 + 0.24 


1.76 + 0.26 


-0.17 


0.16 


0.115 


7.11/7 


1.46 + 0.16 


1.96 + 0.31 


-0.26 


-0.13 


0.121 


4.24/7 


0.82 + 0.11 


0.97 + 0.13 


-0.21 


-0.12 


0.536 


156.00/156 


1.73 + 0.06 


1.77 + 0.04 


-0.20 


0.28 


0.452 


93.70/85 


4.87 + 0.22 


4.87 + 0.22 


-0.03 


0.03 


0.992 


5.03/14 


0.89 + 0.09 


1.00 + 0.07 


-0.20 


-0.14 


0.136 


3.82/8 


1.82 + 0.20 


2.91 + 0.52 


0.46 


0.27 


0.184 


271.57/255 


3.84 + 0.08 


4.39 + 0.08 


0.15 


0.21 


0.914 


43.05/48 


0.48 + 0.02 


1.07 + 0.09 


-0.23 


0.50 


0.124 


155.87/151 


3.85 + 0.05 


20.29 ± 1.18 


0.43 


0.26 


0.722 


31.87/34 


0.49 ± 0.04 


1.35 + 0.08 


0.35 


0.31 


0.665 


54.56/62 


1.87 + 0.30 


23.54+ 14.51 


0.96 


0.02 


0.349 


28.73/36 


0.62 + 0.04 


1.86 + 0.46 


0.58 


0.19 


0.936 


10.79/15 


0.27 + 0.04 


2.85 ± 1.68 


0.68 


0.16 


0.285 


15.82/13 


3.45 + 0.26 


3.99 + 0.33 


-0.15 


0.11 


0.437 


15.50/16 


1.69 + 0.13 
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-0.21 


0.02 
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26.69/20 
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0.14 
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33.88/26 
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-0.08 


0.003 
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-0.02 


0.664 
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0.22 
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0.27 
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0.34 
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0.53 
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0.04 
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3.69/7 
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-0.04 
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0.01 
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0.52 + 0.06 
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0.12 


0.42 
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-0.19 


0.27 
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2.29 ± 1.36 d 


-0.09 
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0.655 
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0.09 
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0.08 
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+0.38 
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+0.08 6 
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TABLE 2 
Continued 
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1 QQ+0.23 


0. 


.09+2-22 

— 0.06 


33.96/27 


1 


.26 ± 0. 


09 


1.53 ± 0.10 


-0. 


20 


0. 


00 





.987 


NGC 


4038 


12 


01 


54 


35 


-18 


52 


10. 


2 


1.18+°?i 

— 0. 19 


0. 


" ' —0.07 


26.56/25 


1 


.98 ± 0. 


15 


2.13 ± 0.13 


-0. 


.02 


0. 


13 





849 


NGC 


4038 


12 


01 


54. 


.77 


-18 


52 


52. 


1 


2.54+0-^ 

0.65 


6. 


■ UJ - — 1.67 


9.15/13 


2 


.93 ± 0. 


69 


17.86 ± 10.75 


0. 


,87 


0. 


07 





517 


NGC 


4038 


12 


01 


54. 


97 


-18 


53 


14. 


9 


1 69+212 

u,y — 0.12 


0. 


21 +0.04 
— 0.04 


84.34/77 


5 


.42 ± 0. 


21 


6.75 ± 0.19 


-0. 


17 


0. 


24 





003 


NGC 


4038 


12 


01 


55. 


65 


-18 


52 


15. 





1.26+°" 

— 0.10 


0. 


03+2-22 

— 0.03 


74.89/84 


7 


.26 ± 0. 


29 


7.64 ± 0.23 


-0. 


,08 


0. 


06 





.807 


NGC 


4038 


12 


01 


56 


.42 


-18 


51 


57. 


8 


1.95+"'" 

— 0.11 


0. 


-,-,+0.03 
— 0.03 


95.67/85 


5 


.09 ± 0. 


19 


6.26 ± 0.20 


-0. 


.21 


0. 


03 





.124 


NGC 


4039 


12 


01 


54. 


.48 


-18 


53 


05. 


.4 


3.23+1-^ 

—0.67 


0. 


19+2?? 

— 0.11 


57.82/14 





.40 ± 0. 


04 


0.85 ± 0.29 


-0. 


.28 


-0 


.19 





.187 


NGC 


4125 


12 


08 


01 


.74 


+65 


10 


34. 


8 


1.51+2f? 

— 0.61 


0. 


.12+2-1? 
— 0.11 


7.50/6 





.84 ± 0. 


20 


0.93 ± 0.17 


-0 


.10 


0. 


13 





.147 


NGC 


4125 


12 


08 


07. 


53 


+65 


10 


28. 


5 


2.01+2-" 

— 0.10 


0. 


Q9+0.00 
u — 0.00 


26.32/38 


5 


.19 ± 0. 


27 


5.88 ± 0.27 


-0. 


20 


-0 


.15 





.208 


NGC 


4125 


12 


08 


17 


.44 


+65 


11 


25. 


3 


-, 21 +0.57 

—0.30 





.03+°" 

—0.03 


7.46/8 


1 


.23 ± 0. 


22 


1.28 ± 0.18 





01 


0. 


.04 





.984 


NGC 


4151 


12 


10 


22 


.37 


+39 


23 


16. 


9 


1.89+0-" 

— 0.18 


0. 


9+0.00 
— 0.00 


4.97/10 


1 


.34 ± 0. 


11 


1.42 ± 0.12 


-0. 


17 


-0 


.21 





112 


NGC 


4258 


12 


18 


43. 


.86 


+47 


17 


31. 


5 


1.81+"-2i 

— 0.28 





1 r+0.10 

— 0.09 


10.71/14 


1 


.06 ± 0. 


08 


1.30 ± 0.11 


-0. 


.22 


0. 


17 





.349 


NGC 


4258 


12 


18 


57. 


.84 


+47 


16 


07. 


1 


2.00+°-?i 

— 0. 19 


0. 


ci +0.10 
u — 0.09 


28.01/32 


2 


.61 ± 0. 


16 


4.22 ± 0.27 


-0 


.14 


0. 


38 





.528 


NGC 


4303 


12 


21 


55 


.46 


+04 


28 


58. 


4 


1.88+2?? 

-0.75 





3o+2-?2 

— 0.19 


6.25/8 


1 


.14 ± 0. 


27 


1.58 ± 0.23 


-0 


29 


0. 


35 





.492 


NGC 


4303 


12 


21 


56 


.57 


+04 


29 


24. 


4 


O 1Q + 1.15 

° °— 0.88 





28+2-2T 

^—0.21 


9.20/6 





.47 ± 0. 


10 


1.11 ± 0.51 d 


-0. 


.43 


0. 


11 


0. 


.147 


NGC 


4303 


12 


21 


58. 


,34 


(11 


■»\ 


12. 


2 


° °— 0.72 





34+2'?i 

— 0. 16 


7.76/10 





.63 ± 0. 


08 


1.71 ± 0.67 


-0. 


.30 


-0 


01 





.748 


NGC 


4303 


12 


22 


01 


.25 


(jl 


29 


37 


4 


^^ — 0.22 


o 


09 +0.00 
UZl — 0.00 


4.41/7 


1 


.64 ± 


18 


1.85 ± 0.17 


-0 


.22 


-0. 


.08 





.016 


NGC 


4365 


12 


24 


26. 


39 


+D7 


16 


53. 


3 


1 oc:+0.44 
0!j — 0.39 


o 


OO+0.21 
< 5 < 5 — 0.17 


4.82/6 


2 


.67 ± 


27 


3.27 ± 0.32 


-0. 


01 


0. 


34 





.289 


NGC 


4365 


12 


24 


28. 


.45 


+07 


19 


05 


9 


2.04+°'^ 

w — 0.35 





09 +0.00 
— 0.00 


2.96/7 





.65 ± 


11 


0.83 ± 0.21 


-0 


.16 


-0. 


.17 





.247 


NGC 


4365 


12 


24 


38. 


29 


+07 


21 


06. 


.8 


3.29+i'S 

^ ^ — 0.94 





14+0.28 

' —0.14 


2.83/5 


o 


.47 + 0. 


10 


0.90 ± 0.54 d 


-0. 


.28 


-0 


.23 





000 


NGC 


4374 


12 


25 


03 


.48 


+ 1 2 


53 


19. 


2 


1.89+0'?? 

—0.52 





o+0.26 
— 0.09 


6.59/8 





.71 ± 


12 


0.83 ± 0.17 


-0 


.28 


0. 


04 





.434 


NGC 


4382 


12 


25 


17. 


17 


+ 1 8 


13 


46. 


.7 


2.24+2^ 
— 0.17 





QO+0.00 

— 0.00 


14.01/16 


2 


.28 ± 0. 


17 


2.68 ± 0.16 


-0. 


20 


-0 


.23 





.413 


NGC 


4382 


12 


25 


20 


.32 


+ 1 8 


13 


01. 


6 


2.96+2^ 

-0.65 





2O+0.23 
° — 0.18 


7.41/10 





.73 ± 


10 


1.59 ± 0.52 


-0. 


.41 


0. 


09 





390 


NGC 


4388 


12 


25 


38 


95 


+ 1 2 


40 


00. 


9 


2.29+°^ 
^ — 0.66 





irr+0.25 
— 0.15 


1.80/5 





.87 + 0. 


17 


1.24 ± 0.32 


-0. 


.25 


0. 


04 





.335 


NGC 


4388 


12 


25 


39 


.10 


+ 1 2 


40 


09. 


5 


2 ^+0.61 
° — 0.50 





.13+2- \% 

^ — 0.13 


9.04/8 


1 


.24 ± 


19 


1.63 ± 0.28 


-0. 


.22 


0. 


02 





.867 


NGC 


4395 


12 


25 


39 


51 


+33 


32 


04. 


3 


^ — 1.26 





45+2 23 

^ — 0.22 


6.52/5 





.07 + 0. 


02 


0.62 ± 0.61 d 


-0 


.31 


-0 


.24 





.100 


NGC 


4449 


12 


28 


17 


.83 


+44 


06 


33. 


9 


1.99+212 

^ — 0.16 





c 7 +0.08 
— 0.07 


54.19/56 





.75 ± 


03 


1.22 ± 0.06 


-0. 


.16 


0. 


41 





.928 


NGC 


4472 


12 


29 


34. 


46 


+07 


58 


51. 


6 


-, 72 +0.24 
— 0.23 





n2 +o.oo 

— 0.00 


2.18/6 


1 


.26 ± 0. 


15 


1.33 ± 0.15 


-0. 


15 


0. 


11 





144 


NGC 


4472 


12 


29 


34 


49 


+08 


00 


32. 


3 


-, o 9 +0.52 

J. . r\ a a 
— 0.44 





22 +0.17 
— 0.14 


4.72/10 





.89 ± 


13 


1.16 ± 0.15 


-0. 


.17 


0. 


13 





.978 


NGC 


4472 


12 


29 


40. 


99 


+07 


57 


44. 


2 


2-7i+sei 

— 0.53 





oy+0.22 
— 0.17 


11.22/9 





.76 ± 0. 


10 


1.55 ± 0.41 


-0. 


,27 


0. 


24 





.260 


NGC 


4472 


12 


29 


42 


.31 


+08 


00 


08. 





1.68+2-52 

— 0.4,5 





-,2+0.15 
— 0.12 


0.24/5 





.93 ± 


13 


1.09 ± 0.12 


-0 


.14 


0. 


01 





.386 


NGC 


4485 


12 


30 


30 


.56 


+41 


41 


42. 


3 


1.08+°-°" 

—0.17 





.16+2-2 4 

w — 0.03 


52.35/54 


5 


.70 ± 0. 


21 


6.34 ± 0.28 


-0 


.13 


0. 


24 





.505 


NGC 


4490 


12 


30 


29 


55 


+41 

-J-<i J. 


39 


27. 


6 


o Qo+0.70 
— 0.58 


2 


rq+0.80 
— 0.62 


21.85/22 





.98 ± 0. 


13 


9.67 ± 5.62 





.10 


0. 


44 





.946 


NGC 


4490 


12 


30 


30 


.82 


+41 


39 


11. 


5 


2.05+°-^ 

— 0.24 


1 


.00+2- \l 

— 0.15 


22.82/28 


2 


.63 ± 


17 


5.10 ± 0.46 


-0 


.01 


0. 


.46 


0. 


.673 


NGC 


4490 


12 


30 


32 


27 


+41 


39 


18. 


1 


1 79+ - 26 

— 0.23 





.40+2, I, 1 , 
— 0.10 


14.71/21 


1 


.69 ± 0. 


13 


2.36 ± 0.16 


-0 


.13 


0. 


37 





.618 


NGC 


4490 


12 


30 


36 


.32 


+41 


38 


37. 


8 


1.93+2-^ 

— 0.24 





,,+0.12 
— 0.10 


19.66/22 


1 


.76 ± 0. 


14 


2.77 ±0.20 


-0 


.11 


0. 


38 





810 


NGC 


4490 


12 


30 


43 


.25 


+41 


38 


18. 


4 


o.8o+2-2« b 

— 0.08 


1 


.04+2-" 
" — 0.15 


35.57/37 


3 


.14 ± 0. 


21 


6.23 ± 0.35 


-0 


.11 


0. 


41 





.761 


NGC 


4486 


12 


30 


47. 


.12 


+ 1 2 


24 


15 


.8 


3.20+2 -52 

— 0.43 





26 +o.n 

w — 0.09 


16.88/18 


1 


.43 ± 0. 


13 


3.25 ± 0.65 


-0 


.32 


-0 


13 





004 


NGC 


4486 


12 


30 


49 


.24 


+ 1 2 
-p j.^ 


23 


34. 


5 


o co+2.84 
— 1.04 





+4+2' T 4 
— 0.14 


17.30/21 





.42 ± 


15 


0.80 ± 1.48 d 


-0 


.19 


-0. 


.44 





.867 


NGC 


4486 


12 


30 


49. 


18 


+ 1 2 


26 


04. 


2 


2.o7+2,- < ;2 

— 0.40 





■06 + 2-hk 

" w — 0.06 


3.47/10 





.90 ± 0. 


14 


1.06 ± 0.14 


-0 


.18 


-0 


18 


0. 


.745 


NGC 


4486 


12 


30 


50 


.79 


+ 1 2 


25 


01. 


.8 


1.94+2-52 

—0.42 





ns +0.12 

. UO r\ r\c, 

— 0.08 


16.67/18 


1 


.26 ± 


15 


1.48 ± 0.15 


-0 


.14 


-0. 


01 





473 


NGC 


4486 


12 


30 


50 


90 


+ 1 2 


23 


25 





4.26+2-Ik 

" — 2.06 





.40+2-S5 

— 0.37 


3.77/9 





.26 ± 0. 


13 


1.28 ± 9.42 d 


-0 


.38 


-0 


.14 





.353 


NGC 


4486 


12 


30 


53. 


.24 


+12 


23 


56 


8 


0.16+2-21 c 

—0.13 





OO+0.38 
■ — 0.61 


22.59/14 





.50 ± 


13 


1.74 ± 0.20 


-0 


.10 


-0. 


.10 





.431 


NGC 


4494 


12 


31 


28. 


58 


+25 


44 


57. 


6 


1.76+2-22 

1 " — 0.27 





m +0.00 
— 00 


5.56/5 


1 


.05 ± 0. 


13 


1.16 ± 0.15 


-0 


.16 


-0 


.16 





.921 


NGC 


4494 


12 


31 


29 


.62 


+25 


46 


21. 


9 


2 5Q+°- 67 
uy — 48 





ns +0.'l4 
U8 -0.08 


2.97/7 


1 


.32 ± 


15 


1.87 + 0.33 


-0 


.23 


-0. 


22 





.491 


NGC 


4552 


12 


35 


41 


22 


+12 


34 


51. 


5 


i fi o+0.'57 
loJ -0.49 





ii+O.Ki 
^J^-O.ll 


12.61/10 


1 


.34 ± 0. 


23 


1.57 + 0.17 


-0 


.16 


0. 


07 





.446 


NGC 


4552 


12 


35 


44 


82 


+12 


33 


40. 


.8 


1 icr + 0.15 6 
L - lo -0.73 





oc+0.35 
• 3o -0.75 


8.60/08 





.57 + 0. 


07 


0.93 + 0.14 


-0 


.08 


-0 


03 


0. 


.947 


NGC 


4552 


12 


35 


45 


.77 


+12 


33 


02. 


4 


2 65+ - 54 





97 +(].13 
z '-0.11 


11.64/10 





.82 ± 0. 


.08 


1.51 + 0.24 


-0 


.34 


0. 


11 





.445 


NGC 


4559 


12 


35 


51 


71 


+27 


56 


04. 


1 


o 1o +0.59c 





n9 +0.02 
uz -0.04 


75.81/67 


9 


.41 ± 0. 


79 


9.83 + 0.38 


-0 


.20 


-0 


,14 





.812 


NGC 


4559 


12 


35 


57 


.80 


+27 


58 


07. 


2 


2.28lg;S 





9Q+0.18 
za -0.14 


10.31/7 


1 


.80 ± 


19 


2.89 + 0.46 


-0 


.28 


0. 


22 





.716 


NGC 


4559 


12 


35 


58. 


57 


+27 


57 


41. 


.8 


i-8o±g:iS 





14 +0.05 


58.71/49 


11 


.07 + 0. 


59 


13.60 ±0.58 


-0 


.19 


0. 


15 





693 


NGC 


4579 


12 


37 


40 


.30 


+11 


47 


27. 


5 


1 89+°- 14 





i - + 0.04 
io -0.04 


65.78/58 


16 


.29 + 0. 


63 


20.21 + 0.75 


-0 


.17 


0. 


04 





526 


NGC 


4579 


12 


37 


43 


20 


+11 


49 


01. 


2 


1 69+ - 30 

J..U»_Q 2g 





n9 +0.00 
uz -0.00 


8.63/8 


1 


.52 ± 0. 


19 


1.65 + 0.20 


-0 


.18 


-0 


05 





734 


NGC 


4579 


12 


37 


53 


.89 


+11 


50 


20. 


1 


o oq+0.58 





u ' -0.07 


8.38/10 


1 


.37 + 


21 


1.60 + 0.25 


-0 


.23 


-0. 


.05 
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TABLE 2 
Continued 



Galaxy 


R.A. 
(h m s) 




Decl. 
(d m s 


) 




Spectral" 
Parameter 


n H 
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X 2 /dof 


L 

("1 n39 

^1U 


obs 

erg s _1 


) 


Lint 

fi n39 — 1^ 
erg s ) 


(H - M) 
color 


(M-S) 
color 


Pks 


NGC 4594 
NGC 4594 


12 
12 


40 
40 


00 
01 


.97 
.10 


-11 
-11 


36 
37 


51, 
21. 


2 
2 


9 (=9+0.47 

5.62+149 


0.1810+1? 
0.36i°;f 7 


8.43/9 
7.55/12 






,91 ±0. 
,66 ± 0. 


08 
13 


1.44 ±0.23 
5.73±4.64 d 


-0.30 
-0.22 


0.05 
-0.53 


0.352 
0.096 


NGC 4621 


12 


42 


02 


.99 


+11 


41 


16. 


3 


1 59+° :7i 


21+ ' 24 
u - z± -0.19 


4.93/5 


1 


.00 ± 0. 


18 


1.24 ±0.17 


-0.03 


0.20 


0.938 


NGC 4631 


12 


41 


55. 


56 


+32 


32 


16, 


9 


1 75+ ' 11 

1 -'°-o.io 


16+ 003 


134.36/105 


2 


.42 ± 


10 


2.94 ±0.08 


-0.18 


0.14 


0.599 


NGC 4631 


12 


42 


11. 


13 


+32 


32 


35. 


.8 


9 7Q+0.38 
z - ' y -0.35 


3 qq+0-63 
J - yy -0.56 


44.86/47 


2 


.20 ±0 


12 


13.81 ± 3.96 


0.57 


0.20 


0.228 


NGC 4649 
NGC 4649 
NGC 4649 


12 
12 
12 


43 
43 
43 


32. 
37 
45. 


08 
.27 
00 


+11 
+11 
+11 


31 
31 
32 


18, 
43, 

33. 


5 
7 
.4 


9 OO+0.76 
z -°°-0.62 
i qe+0.25 
l-»O_0.22 
i co+0.28 
±.oo_ Q 26 


K7+0-38 
u - D '-0 .27 

02+ 00 
u ' uz -0.00 

02+ 00 
u - uz -o.oo 


3.12/6 
4.91/7 
14.13/11 







.86 ± 0. 
.97 + 0. 
.90 ± 


14 
11 
11 


1.73 ±0.47 
1.04 ±0.12 
0.93 ±0.14 


-0.16 
-0.15 
-0.10 


0.40 
-0.26 
-0.10 


0.164 
0.016 
0.986 


NGC 4697 
NGC 4736 


12 
12 


48 
50 


46 
52 


.83 
.72 


-05 
+41 


48 
07 


53. 
19. 


5 



1 68+ ' 37 
1 -°°-0.33 

1 82+ ' 12 


14+° 12 

04+ ' 03 
u - u *_0.03 


14.21/9 
83.76/81 


1 




.00 ± 0. 
.86 ± 


10 
03 


1.16 ±0.11 
0.95 ±0.03 


-0.25 
-0.15 


0.06 
-0.13 


0.246 
0.086 


NGC 4945 


13 


05 


21 


.94 


-49 


28 


26 


6 


1 09+0 06 i> 
± - vv -0 .06 


44+0-04 
u '^*-0.04 


97.11/110 





.77 + 0. 


02 


1.04 ±0.03 


-0.16 


0.41 


0.458 


NGC 4945 


13 


05 


32 


89 


-49 


27 


31. 


1 


i si+o- 12 
1 - 81 -o.n 


76+ 08 
u - ,D -0.07 


113.59/101 





.95 ±0 


03 


1.52 ±0.05 


0.02 


0.39 


0.187 


NGC 5055 


13 


15 


19. 


.55 




03 


02 


3 


9 n +0.15 
z - o± -0.14 


n 10+0. 04 


96.01/83 


13 


.45 ± 0. 


19 


21.28 ±0.95 


-0.29 


0.04 


0.011 


NGC 5055 


13 


15 


39 


33 




01 


53. 


.4 


9 ori+0.58 
^■•5^-0.49 


54+0- 26 
u '°*-0.20 


2.91/5 





.57 + 


07 


1.07 ±0.21 


-0.31 


0.46 


0.780 


NGC 5055 


13 


16 


02 


.27 


+42 


01 


53. 


6 


O CT9+0.37 

°-° -0.32 


56+° 11 


21.85/23 


1 


.29 + 0. 


09 


5.56 ± 1.06 


-0.49 


0.36 


0.000 


NGC 5204 


13 


29 


38. 


61 


+58 


25 


05. 


6 


o qq+0.25 
z.»»_ 22 


10+ 005 


45.10/41 


1 


.20 + 


05 


1.90 ±0.17 


-0.26 


-0.28 


0.443 


NGC 5194 


13 


29 


43 


29 


-1-47 


11 


34. 


.7 


Q 11+0-04 6 
u - ±1 -0.03 


22+° 18 
u - zz -0.30 


3.64/6 





.23 ±0 


10 


1.65 ±0.20 


-0.05 
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a Powcr-law index, T, unless otherwise noted. 
^Disk-blackbody model inner disk temperature, kTi n , in keV . 
c Thermal emission-line model temperature, fcT, in keV. 

^Intrinsic luminosity <10 39 ergs s _1 when fit using a thermal emission-line model. 
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Linear Correlation Coefficients" 
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a () denote marginally significant correlation, — denotes no correlation 



